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Abstract 

We prove the existence of a solution to an equation governing the number density within a 
compact domain of a discrete particle system for a prescribed class of particle interactions taking 
into account the effects of the diffusion and drift of the set of particles. Each particle carries a number 
of internal coordinates which may evolve continuously in time, determined by what we will refer to as 
the internal drift, or discretely via the interaction kernels. Perfectly reflecting boundary conditions 
are imposed on the system and all the processes may be spatially and temporally inhomogeneous. 
We use a relative compactness argument to construct a sequence of measures that converge weakly 
to a solution of the governing equation. Since the proof of existence is a constructive one, it provides 
a stochastic approximation scheme that can be used for the numerical study of molecular dynamics. 
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1. Introduction 



Under fairly general conditions, set out in section 2.2, we prove the existence of at least one measure- 
valued solution to the following equation: 

f(z)P t (dz) = f f(z)P Q (dz)+ f f ±a s (z)Af(z)P s (dz)ds+ f f b s {z) ■ V f(z) P s (dz) ds 
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for t ^ and / a member of a suitable class of test functions obeying the boundary condition df/dn = 0, 
where n is the outward pointing unit normal. This condition corresponds to the imposition of perfectly 
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reflecting boundaries. In equation €. is the configuration space consisting of a physical domain f2 

and a collection of internal degrees of freedom and P t is the measure representing the number density 
of particles in £. The functions at, b t and H t are rates associated with diffusion and the external 
and internal drifts respectively. The measure It describes a source term for the system and K ti ij are 
interaction kernels describing the interplay of particles within the domain. 

This equation has many potential applications for the quantitative description of molecular systems. 
It generalizes Smoluchowski's coagulation equation and simple diffusion models, e.g., (23- Existence 
theorems for coagulation and fragmentation processes have been developed for spatially and temporally 
inhomogeneous conditions in, e.g., EI El El 13 El (see also the references therein). The case of 'cluster 
coagulation', extending this work, where the coagulation rate may have non-trivial dependence on some 
internal coordinates has been examined in [§| . For a survey of some of the known exact solutions and 
a discussion of the mathematical and physical reasons for studying the coagulation equation see the 
review article jTHj . 

Diffusion processes with boundary conditions have been studied using SDE techniques by, e.g., 
[Ill 1121 [T5] . However we do not follow this route since here it is our aim to produce a constructive proof 
for the existence of solutions to equation l|l.l|l . This is in order to provide a Monte Carlo approximation 
scheme capable of simple computer implementation that can simulate the underlying physical system 
accurately. 

The existence problem for combined diffusion, coagulation and fragmentation processes has received 
some attention in |141 1151 1161 ITT) for the case of spatially and temporally homogeneous diffusion co- 
efficients. Van Dongen discusses the case when the diffusion coefficients are non-zero and equal for 
masses no larger than some given mass and zero for larger masses. As a separate case he considers 
the situation when all diffusion coefficients are equal. Collet & Poupaud found a local existence result 
in the former case for a restricted range of coagulation kernels and a global existence result in the 
latter case. Guias, ^JEUi developing the work in [T2|, has presented a numerical approach to this 
problem by dividing the domain into cells then solving each cell as a homogeneous coagulation system 
with appropriate transfer of particles between neighbouring cells; thus combining a finite element and 
Monte Carlo approach. Guia§ provides the necessary mathematical theorems to justify this strategy. 
In contrast to this combined finite element and Monte Carlo method, our resulting numerical approach 
is purely stochastic in nature, and the diffusion, drift and other rate coefficients may have spatial and 
temporal dependence, we also allow for more general forms of particle interaction to take place. 

The paper is organized as follows. In section 2 we give definitions of the physical domain and 
configuration space, and describe the hypotheses we will be making. Then we introduce the class of 
test functions and describe the boundary condition they obey. Before beginning the mathematical 
assault on the problem, we relate the boundary condition on the class of test functions to appropriate 
Neumann boundary conditions on the particle number density when such a density exists with respect to 
Lebesgue measure and certain extra differentiability conditions are assumed. Additionally in section 3 
we characterize the boundary conditions associated with the internal degrees of freedom. 

Section 4 is devoted primarily to the proofs of proposition 14.151 and theorem 14. 1 71 that are needed 
in section 5. Proposition 14 . 1 51 provides a bound for an appropriate stochastic generator associated with 
the diffusion and drift terms in equation (|l.l(l . whilst theorem 14 . 1 71 provides a result that allows us to 
deduce the uniform convergence of a discrete approximation to the derivative terms in the integrand of 
the RHS of equation (|l.lfl . To implement the diffusion process each particle's position is displaced by 
a multivariate normal random variable that is scaled proportional to the square root of the time step. 
The resulting displacement is modified to account for the drift terms and for any reflections off the 
boundaries. Since in a finite time the distance travelled by such a particle becomes unbounded as the 
time steps decrease, the proof relies on the cancellation of motions in opposite directions and it is this 
that requires delicate analysis, particularly in a neighbourhood of the boundary. We point out that the 
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use of such a multivariate normal random variable makes direct contact with the mathematical definition 
of Brownian motion and with the paper of Roth |20j on generating a Feller semi-group associated with 
a simple diffusive process. 

In section 5 we formulate a relative compactness argument for a sequence of stochastic processes 
taking values in an appropriate product space. The components of each term in the sequence comprise a 
measure representing an approximation to the particle number density at each time and a non-negative 
real number approximating that time. This fictitious time is expedient for recasting equation HI. If) as a 
pair of equations with no explicit time-dependence. Tightness of the sequence is proved and Prohorov's 
theorem is used to deduce the existence of a convergent subsequence. We verify that the limit is a 
solution to equation in the weak sense when the Linear Growth Condition holds and in the vague 
sense if the growth condition on the internal drift and self-interaction kernel are relaxed. Finally a new 
sequence of processes is presented with the same jump chain of measures but with modified holding 
times determined by the the fictitious time variable. This has the effect of removing the fictitious time 
quantity from the construction. This final result, theorem l5.31l provides a simple approximation scheme 
that is in essence a simple extension of that of Marcus [21] and Lushnikov [22| except that the holding 
times are determined deterministically from the current state. We remark that proposition 15.291 and 
theorem 15.311 can be modified simply to remove all spatial dependence and in this way we recover the 
earlier results of UHSj. 

In section 6 we provide an explicit description of the approximation scheme for generating sample 
paths for the stochastic processes, while in section 7 we make some remarks about potential applications 
to physical situations. These include a brief discussion of the implications and limitations of the various 
hypotheses assumed in the proof. Appendices A, B and C contain important supporting results that 
do not fit naturally within the main text. 



2. Definitions and Hypotheses 

For any pair of Banach spaces A and Y, let Cb(X; Y) denote the bounded continuous functions from 
X to Y and write Cb(X) = C{,(A;R). Let C n (X;Y) denote the n-times continuously differentiable Y- 
valued functions on X. We denote by C C (A) the compactly supported real-valued continuous functions 
on X and define C™(X) — C n (X) H C C (X). Additionally Mb{X) denotes the space of bounded non- 
negative Borel measures on A, 23(A) is the Borel cr-algebra of X and M(X) is the space of bounded Borel 
measurable functions on X. For any P 6 A^s(A) and / S M(X) we will write P(f) = j x f(x) P(dx). 



2.1. The Configuration Space € 

Definition 2.1 Let fl = o; _1 (— oo,0) be a bounded open subset of the Euclidean inner product space 
(M. 1 , ■), with the induced norm denoted by || • ||, and where lo € C 2 (R dl ) with ||Vw(a;)|| > 1 for all 
x G dfl = cj _1 {0}. Let Q be the closure of ft and write V = (pi, . . . , d^) for the gradient operator on 
R dl and A for the Laplacian. 

Definition 2.2 Define the expressions 




and 




(2.1) 
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where co(A) denotes the convex hull of the set A. 

Definition 2.3 Let a™^ ^ f>^" for m G N and /i = 1 , . . . , c?2 • Define the sequence of compact sets 
r m = 11^=1 [ a p™^ ' ^ m ^] R^ 2 • W^e endow T m with the Euclidean inner product induced from R d2 and 
write a dot for this inner product and \\ ■ || for the Euclidean norm. In addition we let V = (di, . . . , dd 2 ) 
denote the gradient operator on M. d2 . 

Definition 2.4 Define the configuration space £, and its extension both subsets o/N 
by 

OO u OO 

£ = U ({m} x n x T m ) and £ # = U ({m} x co(ft) x T m ), (2.2) 

m— 1 m—1 

and endow each with the Euclidean metric induced from R 1 + d i+" 2 . Define the three restrictions: n m : 
€ -> N, tt^ : £ -> O and 7r r : £ -> M d2 6y 

^(zo^ij.-.j^+dj = z ; (2.3) 
7i^-(^o,«l)---)«d 1 +d a ) = (zi,---,z dl ); (2.4) 
7r r (z ,zi,...,z dl+d2 ) = (z dl+ i, . . . , z dl+d2 ), (2.5) 

so that 7Tr(z) € r^ m ( a ). It is useful to define eo, ei, . . . , e^, ei, . . . , ed 2 to be the canonical basis for 
jga+di+d 2 j n ^ e o\) V i ous -way. For convenience we will refer to n m (z) as the mass of a particle with 
coordinates z. 



2.2. Hypotheses 

In this subsection we present two choices of physical hypotheses that are sufficient to prove an 
appropriate convergence result. 

Hypothesis 2.5 Suppose the source rate is given by the continuous map I : [0, oo) — ► A4b(£), with 
1 : 1 1— > I(t) = I t satisfying the boundedness condition 

A( IV2 > = sup {/ t (^ lv2) ) : < t < oo} < oo (2.6) 

where we define tt^P (z) — (ir m (z)) q and where X some the non-negative integer defined implicitly by 
equation 1)1.1(1 . 

Hypothesis 2.6 Suppose the diffusion coefficient a : [0, oo) x £ — > K (with a(u,z) = a u (z) = o- u (z) 2 ) 
is continuous and there exists a bound o~oo such that 

< a u (z) ^ a x (2.7) 

for all u £ [0, oo) and z£C. Write 

cr^ = 1 A inf{er t (z) :z££, 7r m (z) < to, < * < T}, (2.8) 
which is strictly greater than zero by continuity on the compact set [0, T] x 7T" 1 !!, . . . , to}. 
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Furthermore assume that the drift function b : [0, oo) xC-* R dl (with b(u, z) = b u (z) ) is continuous 
and obeys the boundedness condition ||o„(z)|| ^ for some b^ € M and additionally satisfies the 
boundary condition 

b u (z)-n(n Ti (z))=0 where n(x) = (2.9) 
/or a/Z it S [0, 00) and z € £ swc/i f/iai 7Tq(z) S 9f2. 

Hypothesis 2.7 Suppose the internal drift rate H : (t,z) 1— > H t (z) = (H tt i(z), . . . , H t ,d 2 { z )) is a 
continuous function from [0, 00) x £ to M d2 where 

H™ = 1 Vsup{||i/ t (z)||/7r m (z) : z S £, 7r m (z) m, ie[0,oo)}<oo (2.10) 

/or oZi m£N, ^4s a matter of notational convenience we set = 1. We will also impose the boundary 
conditions: H ttfl (z) whenever (71^(2))^ = aj^ m ^ z ^ and H t ^{z) ^ whenever (tt-t(z))^ = oj? m . 

The purpose of these boundary conditions is to ensure that the internal drift does not take particles 
outside of the domain of the internal coordinates. 

Hypothesis 2.8 Suppose there exist uniformly continuous interaction kernels K% t j for i ~ I,..., I, 
j = 1, . . . , J and (i,j) 7^ (1, 1) with K lyj : [0, 00) x <T — * Mb{&) — ®k=i -^s(^) where we write 

K u>i ,j(zi,...,Zi,Bi x---xBj) = K itj (u,zi,...,z i )(B 1 x---xBj) (2.11) 

which we require to be totally symmetric in Bi, . . . , Bj 6 23(£), totally symmetric in Zi, . . . , zi and to 
satisfy 

i j 

y^^m(zi) = y^^m(w n ), K u>iij (zi, . . . , Zi, du?i, . . . , dwj)-a.e. (2-12) 

1=1 n=l 

for all Z\, . . . , Zi £ £, i = I, ... ,2, and j = 1, . . . , J , i.e., the interaction is mass preserving. Moreover 
we assume that there exists a constant Kao such that 

-^—K Uy i. J (z l7 . . . , Zi,<L 3 ) < ifoo7rm(zi) • • -iT m (zi) for ^ (1, 1). (2-13) 

In addition we assume that there exists a continuous ( though not necessarily uniformly continuous ) self- 
interaction kernel, Ki t \ : [0, 00) x £ — > A4b(£) where we write K Ut i t i(z, B) = K\ t i(u, z){B) satisfying 
ftm{z) — TT m (w), K Ut i_i{z,&w)-a.e. and possessing the bound: 

K™ = IV suv{K txl (z, £)/7r m (z) ir m (z) ^ m, t S [0, 00)} < 00. (2.14) 

Analogous to our definition of we set = 1. 

Hypothesis 2.9 We assume that the initial measure-valued particle density has finite mass: 
00. 

Remark 2.10 If both and H!^ are bounded functions on N then we say the system obeys the 
Linear Growth Condition. In which case we may chose and Hoo such that for all m £ No we have 
K"^ < Koo, < iJoo and equation (|2.13|l still holds. When this additional condition holds we will be 
able to prove that equation (f 1 . 1 1> holds in the weak sense. When it does not we show that equation (|l.l(l 
is true in the vague sense. 
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2.3. The Class of Test Functions 

Definition 2.11 Define the normed space of junctions 

3 = {/: /(m,0 eC 2 (Hxr m ), for all men, |||/||| < oo 

and Vw(7r ?T (z)) • V/(z) = for all z G € such that K^(z) G 90} (2.15) 

where 3 and its subspaces have the norm ||| • ||| defined by 

HI/HI = max{||/|U, HV/lloo, IIV/IU, HW/lloo} (2.16) 

IIV/IU = sup ( £ |a,/(z)| 2 j ; ||VV/|U=sup( 53 IW(*)| 2 | (2.17) 
\i=i / \i,i=i / 

and 

IIV/IU = sup El^/WN ' ( 2 - 18 ) 
Furthermore we define 3c — 3 n C c ((£) and Zei £/ie map i? : 3" — ► No fee given by 

= fmax{7r m (supp/)} iffedc; ^ 
1 otherwise. 

Remark 2.12 By proposition IA.1I (3c, || • |||) is separable, therefore there exists a dense sequence 
fn S 3c- We set /o = 1. Proposition IB . 21 implies that for each / € 3c U{/o} there exists a second-order 
continuously differentiable extension /# : — > K with the property that |||/^|||* is finite, where we 
define the mapping ||| • |||# on each extension by equations H2.16tl2~TB|) with €. replaced by £"\ 

The reason for introducing these extensions is to allow us to assume that the test function is defined 
on a convex domain for each integral mass. This allows us to make elementary estimates using Taylor's 
theorem applied to / # without regard for the possibility that the straight line segment between the two 
points of interest may not lie entirely within the domain of the unextended test function /. In this way 
we need not assume any convexity property of the physical domain Q. 

Definition 2.13 Let the sequence (f n )^Li € 3c be a dense subset of (3c, ||| • |||) with the following 
properties 

|||/# |||# < 2"; H% fn) K% fn) ^2i» (2.20) 
and (f n )^—i is closed under the operation f t— ► f(')H-ir m (.)<,k f or a ^ G N- 

Remark 2.14 It follows from definition 12 . 131 that for every / € 3 there exists a subsequence of /„ such 
that 

bp-lim f nk = /; bp-lim 9j/ nt! = (%/; bp-limd,dj/ nfc = d t d 3 f (2.21) 

k — >oo k — >CO k — too 

for i, j = 1, . . . , di and 



bp-lim 9 M /„ fc =d„f (2.22) 
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for n = 1, . . . , g?2 (where bp-lim denotes bounded pointwise convergence, i.e., bp-lim/„ = / if and only 

n — >oc 

if f n {z) — > /(#) as 71 — ► oo for all z S £ and supsup |/n(z)l < oo). Moreover if / S 5c then we may 

n ze<t 

assume that 

R(fn k )^R(f) (2.23) 

for all k. 

Definition 2.15 We endow A4b{£) and its subspaces with the complete metric generating the weak 
topology given by 

oo 

dwcak(P,Q) = 5>-"(|P(/„) - Q(fn)\ A 1) (2.24) 
and where fo = 1. FFe wiie P„ — P to signify weak convergence. 

Remark 2.16 We note that since £ is separable Mb^) — ®j = i-^s(£) an d if — 1 P then 

Finally, to establish our conventions, let us make note of the following definition and the subsequent 
remark; they will play a frequent role in proving the convergence of various quantities in the rest of the 
paper. 

Definition 2.17 If (X, dx) and (Y, dy) are metric spaces and f : X — > Y is a continuous function we 
define the modulus of continuity by 

w(f,S)= sup d Y (f(x),f(y)) (2.25) 

d x (x,y)^S 

whenever the RHS exists. 

Remark 2.18 The function / is uniformly continuous if and only if 

limu;(/,<5) = 0. (2.26) 

In particular if / : X — > R m is continuous and compactly supported then / is uniformly continuous and 
equation l|2.26|l holds. 



3. Neumann Boundary Conditions on the Number Density 

In this section we relate the properties of weak/vague solutions to equation (|1 . 1|> to the boundary 
conditions on the number density with respect to Lebesgue measure when such a quantity exists. Let 
us write z = (m,X) and X = (x,X) with m £ N, x £ Q, X S T m and assume that the following 
definitions can be made with respect to Lebesgue measure: It({m} x dX) = Jj™(<Y) dX, 

K tjid ( Zl , ..., Zi , {„:} - .;>•; K,| x dy,-) = K ™}<~w< n i>-< n i :.V; x h y, . . . . , y 5 ) dy! . . . dy,, 

(3-1) 

where K t ' " '" 3 and I™ are continuous functions on their respective domains, P t ({m} x dX) = 
c\ n (X) dx\ with a t S C b 2 (£), b t S C^C; R dl ), H t e C£(£; R d2 ). Suppose too that cf (X) is continuously 
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diffcrcntiable in time and twice continuously differentiable on f2 x r m , then we discover for all X G 
x (r m )° (the interior of f2 x T m ) that 

-^1 = ^A(a t (z)cT(X)) - V • (b t (z)cT(X)) - V ■ {H t {z)cf{X)) + I?(X) 

+ EE fi77~ TTi £ £ d ^"-£ dy^ f_cr(x 1 )dx 1 --- [c^ix^dx, 

Z — 1 7 — 1 t "1 j — l 1 H 

J mi , . . . , m — 1 

ni ,. . . — 1 

X K ^ r ^u...,n J . 1 ,m (x ^ V( J| ... ..J, ..V, 

,_i ^"J-J-j- JaxL, ./nxiv ^nxr mi Jnxr m . , 



2=1 7 = 1 

J ni ,. . .,nj — 1 



x JQ.-™^^, V, :..V. J';... (3-2) 

This follows by considering a suitable sequence of test functions concentrated at z and compactly 
supported in £1 x (T m ) . Now multiply equation (|3.2fl by f(z) £ 5 C , subtract the time derivative of 
equation and apply the divergence theorem to leave the boundary terms: 

= / / [y(z)\7(a t ( Z )cT(X))-^a t (z)cT(X)\7f(z)-f(z)cT(X)b t (z)]-n(x)dSdX 

"E l,^"' l<7 X ^ fTM^--- t<\ AXd * I [HtA^T(X)f(z)] X /Z bl C dx(3.3) 

~,Ja\ m) Ja (m \ Ja {r l\ J a ( , m) 

where n(x) is the outward pointing unit normal at x. The last two terms in the first integrand on 
the RHS vanish by the definition of the boundary condition on the test functions and by hypothe- 
sis 12.61 respectively. By taking a suitable sequence of test functions concentrated at z we deduce that 
d{a t (z)cl n (X))/dn = for x e dQ and c™{X) = whenever = a^ n) or b { ™ ] and H t ^(z) ^ 0. 



4. Perfectly Reflecting Boundaries 

The diffusion, internal and external drifts of each individual particle are to be approximated by 
a sequence of discrete jumps in the particle's position in the configuration space £. To complicate 
matters these jumps may involve reflections off the boundaries of the physical domain. In this section 
we provide some detailed calculations giving estimates and convergence results for our diffusion and 
drift approximations which will be vital in section 5. 

Proposition 4.1 (a) If x,y € dfl with \\y — x\\ =9^0, and n(x) = — |— jjy then there exists a 
constant Bq such that 

\n(x) ■ k\ < B 6 where k = . V ~ X . (4.1) 

\\y-A\ 

(b) Ifx,y £ fl with \\\?uj(x)\\ ^ 1 and Vw(j/) ^ then there exists a constant Aq such that 

\\n(y) - n(x)\\ ^ A 9. (4.2) 
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Proof, (a) Use Taylor's theorem with y = x + dk and the facts that uj(x) — tu(y) — and ||Vu;(x)|| > 1 
to deduce that 



\n(x)-k\ ^ \9^ 1 uj(x + 9k)-e- 1 Lu{x)-VLu{x)-k\^Bo<: 



where B = |||VVw|| c 
(b) Compute 



\\n(y) - n(x)\\ 



Vw(y) Vw(i) 



||Vw(y)|| ||Va,(x)| 
\\Vu>(x)\\ - ||Va,(t/)|| 



Vw(y) 



Vcj(y) - Vw(i) 



||V W (i)|| ||Vw(»)|| w \\Vcj(x)\\ 
||Vw(a:)||-||V«(tf)||| + ||Va;(»)-Vw(a:)|| 



||Va,(x)|| 
< 2\\Vui(y) - Vuj(x)\\ ^ A 9 

where A = 2||Wo;|| 00 . 

Definition 4.2 Let £ : x R dl ft &e de/med fry 

£(x, /c) = a;*, for some x* such that \\x + k — x*|| = min{||:r. + k — y| SI}. 
T/iis exists by the compactness o/f2. 

Proposition 4.3 If x £ dSl, k ■ n(x) = and \\k\\ < \B^ X then 

\\i(x,k)-x-k\\ < 2 J B ||fc|| 2 . 

Proof. Consider y = x + K where K = k — 2£?o||Vaj(x)||~ 1 ||fc|| 2 n(x), then 

u{y) ^ (k-2B \\VLu{x)\\- 1 \\k\\ 2 n(x))-Vuj{x) + ±\\VVLu\\ oo \\K\\ 2 
< -Bollfcll^l-^Ufcll 2 ) <0, 

and thus y £ £1, it follows that 

\\£(x,k)-x-k\\ ^ \\y-x-k\\ = 2B ||Vw(x)||" 1 ||fc|| 2 ^2B Q \\k\\ 2 

as required. 



(4.3) 



(4.4) 

(4.5) 

(4.6) 
(4.7) 
I 

(4.8) 



(4.9) 



(4.10) 
(4.11) 



(4.12) 
I 



Definition 4.4 Define the map 7 : fl x 



O by the following procedure. If x € Q, and k £ R 1 \ {0} 



then let xq = x and ko — k/\\k\\, let i9„ ~ inf{i € [0, ||fc|| — J2" =1 ■ x " +tk n £ £1}, x n +% = x n + r2 n fc n , 
kn+i — k n — 2(k n ■ n(x n+ i))n(x n+ i) . We set j'(x,k) = linin^oo x n (this exists since the {x n } are 
contained in the compact ball of radius \\k\\ centred at x, so there exists an accumulation point, which 
is unique since \\ x n ~ ^n-lll = Y^=o^n is finite). We define 7'(x,0) = x. Finally we set 

j(x, k) = "f'(x, k) if Ys^Lo^n = \\k\\ otherwise let j(x,k) = ^(-/'(x,k),(\\k\\ - Y,T=0^r)k oo ), where we 
remark that proposition ^. 1\ implies that for n > m, \\k n — k m \\ ^ 2^™ =m+1 | fc^ • n(a;i ) | 2Bq 5^™_ m+1 
which shows that the k n are Cauchy and the completeness of the (di — 1)- dimensional Euclidean sphere 
implies that k n — » koo for some fcoo with ||fcoo|| = 1 and 71(7' (x, k)) ■ fcoc = 0. 
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Definition 4.5 Define for r > 0, 



dQ r = {yGfl:d{y,dn) <r}. 



(4.13) 



where d(y, A) = inf{||y — <z|| : a € A}. 



Proposition 4.6 There exists < S < \B Q 1 such that for all x,y € dQ 2S with \\x — y\\ < 25 we have 
\\n(x) — n(y)\\ ^ ^oll^ — y\\ and if additionally x,y £ dfl then \n(x) ■ (x — y)\ ^ -Boll^ — y\\ 2 ■ 

Proof. Since ||Vw|| > 1 whenever x <E dil, d£l is compact and x i— ► ||Vu;(x)|| is continuous, it follows 
that there exists r] > such that ||Vw(x)|| ^ 1 + rj for all x 6 d£l. The uniform continuity of Vw 
on il implies that there exists So > such that if ||x — y\\ < So then ||Vw(x) — Vto(y)\\ < r] and 
thus ||Vw(y)|| ^ ||Vw(a;)|| - ||Vw(a;) - Vw(y)|| > 1. Let 25 = S A B^ 1 and the result follows by 
proposition 14. II 1 



Proposition 4.7 For all z € C with n^(z) — x, k such that d r = \\k\\ and j(x,k) ^ X2 we have 
that 

|( 7 (a;, k) - x - k) ■ Vf(z)\ < 2(2£? + A )\\ VV/ lUlfcH 3 (4.14) 

/or eac/i / £ 5- 



Proof. Suppose that Xi, . . . , xjq £ dtt with 7(2;, k) — xjv+i and iV > 1 then 



Mx,k)-x-k)-Vf(z)\ = 



N 



^^ r (fc r -fc )-V/(; 

r=l 
N r 



(4.15) 
(4.16) 
(4.17) 



r— 1 p— 1 

Af r 

= 2 EE^(v »(*p))n(»p) • [ww - v /m 

r— 1 p— 1 

N r 

< 2 J B ||VV/|| oo ||A;||^^^ p + 2||VV/|| oo ||A : ||^|fc j vn(a; J v)| (4.18) 

r— 1 p— 1 

< 2B ||VV/|| oo ||£ : || 3 + 2|| VV/IUIIfcll^jvdfcjv-i ■ n(xjv-i)| + >Mjv-i) 

(4.19) 

s=: 2(2B + J 4o)||VV/|| 00 ||fc|| 3 (4.20) 
where z p = (-7r m (z), x p , 7rr(z)). We have made use of proposition ^, ll to derive inequality (|4.19(l and that 

||V/(z) - Vf(z p )\\ sC ]T ||V/(^_!) - V/(z r )|| < ||VV/||ooX> < HVV/Hoollfcll, (4.21) 

r— 1 r— 

noting that the straight line segment [z r _i,z r ] is contained within £ to derive inequality 14.18|) . In 
the event that no such N exists the above proof shows that |(xjv+i — x — k) ■ V f(z)\ ^ 2(2Bq + 
Ao)\\ VV/Hoollfcll 3 for all N > 1 and the conclusion follows from 7(2;, k) = "f'(x, k) = limAr^oo xn- 1 
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Proposition 4.8 (a) If S is as in proposition \4-b\ x £ dVl s and \k ■ n(x)\ > (Aq + i?o)||fc|| 2 then 
j(x,k) = x 2 . 

(b) Ifj{x, k) = x 2 and x = tTq(z) for some z£C then \(^(x,k) - x - k) ■ V/(z)| ^ 2|| VV/HooUfcH 2 for 
each /ej. 

Proof, (a) If x\ dtt then 7(2;, k) — x + k = x\ — x 2 and there is nothing to prove, so we assume that 
x\ 6 dil. It follows by proposition 14. II that || 77.(3:) — ^ ^-oll^H and we have 

%i = k -2(k -n(x 1 ))n(xi), (4.22) 

so that 

\ki ■ n( Xl )\ = \ko ■ n( Xl )\ ^ \ko ■ n{x)\ - \\n{x) - n{x x )\\ > B \\k\\. (4.23) 

Proposition ^. II also implies that 7(2;, k) = 7(0:1, (||fc|| — $0)^1) — x\ + (\\k\\ — $o)ki — x 2 ^ dQ. 
(Jo) We have 

\{-y(x,k)-x-k}-Vf(z)\ = |{(||*||-i?o)(fci-^)}-V/(z)| (4.24) 

= \2d 1 {k -n{x 1 ))n{x 1 )-{Vf{z)-Vf{n m {z),x 1 ,TT T {z))}\ (4.25) 
«S 2||VV/|| 0O ||fc|| 2 (4.26) 

where we have made use of the boundary condition n(x\) ■ V/(7T m (z), x\ 1 tty(z)) = and that the line 
segment lies within f2. 1 

Definition 4.9 Define the map jj^ : [0,oo) x £ — > R d2 associated with the internal drift to be 

(Jm(«, *))„ = 4" m(z)) V M*)/« + c N 2 HuAz)) A ^'" (z)) for p = 1, . . . , d 2 . (4.27) 

Proposition 4.10 If f £ $ c , — > s and cjv — > 00 as N — > 00 i/ien 

firo(«f, *) ~ M*) ~ ^) ' V/(z) - (4.28) 

uniformly on £ as N —f 00. 

Proof. We may assume that 7r m (z) ^ -R(/) as otherwise the terms in the sequence are identically zero. 
By the boundary conditions on the internal drift ( hypothesis 12. 7|) and equation (|4.27(l . if c%(j^,(u^ , z) — 

7r r (*)V # then either Jf>f , z) ^ <? N (jmW ,*) " and (j£(uf ,z)) M = 6^- W) 

letting us deduce that iZ^uf ,z + (j^,(u^ , z) - n T (z)) ^) ^ 0, or ff M (nf , 2) < Cff(jm( u ^ > z )-^r{z))n 
and (jm^f ,z)) M = a|T m(z)) leading to H^{u^ , z+{j^(u^ , z)-irr(z))^e ft ) ^ 0. Let <? M (A) = H„{v%,z+ 
A(j^)(u^, 2) — 7^(2))/^) — c 2 N (j^(u^ , z) — 7Tr(^))p then : [0, 1] — > K is continuous and additionally 
3,11(0) and 3^(1) are either zero or have opposite sign. Therefore there exists A M S [0,1] such that 
<^(A M ) = 0. Take A„ = 1 if c 2 N (j"( u »,z) - 7* (*)),* - H^,z)._ Set C (ai) - a *) - 

7rr(«)) M e M ), then we have that - z\\ sC \\H(u*f , z)l,r m llc^ 2 < H™ (f) R(f)c N 2 and 

V C N / 
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< \\Vf\\ooJ2\ c Nti&(<,z)-Mz))»-H»(u?,z)A WmizHR(f ^ (4.29) 

M = l 

< HV/IU^I^^f,^)!^^))^/) -^(«^,«)l ffmW <ii(/)| (4-30) 

^ ||V/|| cc ^ W (i7 M (-,-)a. m (.)^fl(/),kf -5| +J ff- (/)J R(/) Cw 2 ) (4.31) 

which tends to zero uniformly on £ as AT — > oo using the uniform continuity of H^{-, ')^-ir m (-)<ii(/) on 

[0, s + 1] x £ and the definition of the modulus of continuity (equation l|2.25[l ). 1 



Proposition 4.11 If Z ~ N(0,I dl ), n S R dl roii/i ||ra|| = 1, Z ^ 0, m > 1 and cat —> oo as TV 
i/ien 1/iere exist constants D{ m and -D 2 ' m smc/i 1/ial 



oo 



^ Z (\\n\ z . nKc ^ (A+Bmm) ) <£ 1 -T_ 1 a . (4.32) 



Proof. If rfi > 1 then write ||Zj| = 7' and Z ■ n — rsin#, then compute: 
E (ll^ll a |z.n|<cJ v r m (A+B||Z||"')) 

2 2-id! /-oo / .sin- 1 ((c^ m (Ar- 1 +Sr m - 1 ))Al) 

= — — j =— ■ / / cos rfl - 2 6id6i r '+«ii-i e - 3 r dr (4 33) 

V7rr(±di - 5) Jo Jo 

r- ^— ^ HiAr- 1 + Br^y+^e-^ 2 dr (4.34) 

fr- l T$di - |)c™- x Jo 

z_i - 1)) + 2^gr( 



2* 



2W0F(,4r(l(Z + dl - l)) + 2^BT{\{1 + m + d, - 1))) 

^ wTj H _m-l ' (4.CS0J 



whereas if dj = 1 and I > then 

~2 r °° 



^ Z (\\Z\\ l l\z.n\^ N -(A + B\\z\M) = r l e-? r \ rsiA+Brm dr (4.36) 

< — -r\ - / r'e-^'Mr" 1 +Sr m - 1 )dr (4.37) 

V 7T 7 

2*~* (T(|0A + 2?r(% (Z + m))B) 



< V2 (4.38) 



The situation is more delicate when d\ — 1 and 1 = 0. First consider B ^ 0, in this case we use that 

which follows by noting that the function f(r) — c]^ m (^4 + Br m ) — r has exactly two positive roots 
with f{2Ac 1 N m )J(^B-^c N ) ^ provided that 

i l i 

c N > 2 m - 1 A""B m ( m_1 ) . (4.40) 
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To prove this, the inequality 2(1 — 2 ) ^ 1 for X > is useful. It follows that 



^ Z {\\Z\\ l hz-n\^ N - { A +B \ m ) < e-4-dr 



2Aci 



1„2 



- / e-^ 2 (2 J B^Tc w 1 r) m - 1 dr (4.41) 
ttJi CnB - — 

2iA + 2 H ^r(im)S , s 

*S ^_ m _/ 2 ■ (4.42) 



When 5 = the computation proceeds as above except that we note that there is only one zero of the 
function /. The finite number of times when inequality H4.40[l may be violated can be dealt with by 
taking either D l { m or Dlf"" suitably large. I 

Corollary 4.12 If cn — * oo as N — > oo then there exist a function Fi(A) such that 

E Z l| fe .„(*)|<A W < J^ {A \ +1 (4-43) 

°u,7r m (z) C JV 

for all zeC, where I G N , x = fe(z), Z ~ N(0,I dl ) and k = tJu ^ Z + Mil. 

Proof. We make use of a special case of Holder's inequality: \X + Y\ ^ 2'~ 1 (|X|' + \Y\ l ) which follows 
for I ^ 1 by considering the vectors (X, Y) and (1, 1) together with the pair of conjugate indices I and 
y^t. By inspection this inequality also holds for 1 = 0. If \k ■ n(x)\ < A\\k\\ 2 then 

\Z.n(x)\^ a + mr (4.44) 
cn 



where 

+ 2 -4 & oo SU P C JV .> . 2 

a= ■. - and /3 = go0 . (4.45) 

^min ' _uun v 7 

lt,7T m (2:) U,7T m (2) 

The result is completed by the use of proposition 14. Ill I 
Proposition 4.13 The following is a (Taylor's theorem) identity 

F(x + h,y + k)- F(x,y) = [ hD 1 F{x,y + tk)dt+ [ kD 2 F{x + th,y + tk) dt 

Jo Jo 



[ f h 2 D 2 1 F(x + th,y + sk)dtds 
Jo Jo 



(4.46) 



where F : R 2 — > K is continuously differentiate, has continuous second order partial derivative in the 
first argument and the domain of f contains the closed cuboid consisting of the points with coordinates 
(x + th, y + sk) for s, t G [0, 1] . 

Proof. Write t — sr in the last integral on the RHS and perform the integrals. 1 
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Remark 4.14 Proposition 14 . 1 31 implies that if k = — h —5-^, x = ttq{z) and 

Cat c n 

G N (u,z,Z) = f fr(x,k)-x-k]-Vf(K m (z),x,(l-t)Mz)+tjm(u,z))dt 

+ f k-(Vf(n m (z),x,(l-t)ir r (z)+tj&(u,z))-Vf(z)) dt 
Jo 

tJ=1 Jo Jo 

x (d l d ] f*{-K m {z), (1 - t)x + tj(x, k), (1 - s)TTr(z) + sj^,(u, z)) - d l d J f{z)) dtds 

di 

+ ^2 \{"f{ x > k) - x - k)i(^(x, k) — x + k)jdidjf(z) 

+ J \j&(u, z) - 7rr(z)] ■ (v/ # (7r m (z), (1 - t)x + tj(x, k), (1 - t)irr(z) + tj&(u, z)) 

-V/(z)) dt 

+ z) - Mz) - • Vf(z) + £ ^d 3 f(z) (4.47) 

then if Z ~ N(0,I dx ), 

E z (c 2 N G N (u,z,Z)) = E z (c 2 N f(n m (z), 1 (x,k) 7 j^(u,z))-f(z)) 



\a u {z)Af(z) - b u (z) • V/(z) - H u {z) ■ V/(z). (4.48) 



2 



Proposition 4.15 // / G J C U {/o}; c w ~~ * 00 as N — > 00 i/iew i/iere exists some constant L such that 
for all z € £ 

|4E Z (/(^mW,7(^),fc),JiD(«^))-/W)| ^L7r m (2)i/-(/)lll/ # lll # (4-49) 
cr u {z)Z 6„(z) 



/or a/Z AT where Z ~ AT(0, ) and 



Cat 



Proof. If / = /o there is nothing to prove, so we assume that / e 3c. If cc = ttj^z) € <9f2 5 then 
propositions 14.71 and I4.8f b) imply that 

\c 2 N E z ({^(x,k)-x-k}-Vf(z))\ 

< 2c^||V/|| 00 E z (l || fc || >4 ||fc||) +E z ( CAr { 7 (x,/c) - a; - k} ■ V/(*)l 7 (x,*)=* a ) 

+ E z (c 2 N { 7 (x, k) —x — k) ■ Vf(z)t l{x . k) ^) (4.50) 

< 2(<T 1 ||V/|| 00 + HW/lloo + (2S + A ) ( 5- 1 ||VV/|| oo )E z (c Ar ||fc|| 2 ) (4.51) 

< 4(<T 1 ||V/|| 0o + HVV/lloo + (2S + A )(5- 1 ||VV/|| oo )( ( T^d 1 + supc^ 2 ) (4.52) 

whereas if x £ dfl s then 

\c 2 N E z ({ 1 (x,k)-x-k}-Vf(z))\ < 2c Ar ||V/|| 00 E z (lL|| fc || >A -||fc||) (4.53) 

< 2||V/|| 0O 5- 1 E z (c Ar ||fc|| 2 ) (4.54) 

'AT 



< 4||V/|| oo r 1 (a 2 d 1 +fe 2 c supc A ; 2 ). (4.55) 



JV 
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Note that 

\c 2 N E z (k- V/(z))| = \c N E z (a u (z)Z ■ V/(z)) + & u (z) -V/(z)| (4.56) 

= \b u (z)-Vf(z)\ (4.57) 

< ftoollV/IU (4.58) 

so that if we set 

L = 4(1 + J" 1 + (2B + A )<5- 1 )(a^d 1 + 6^ supc w 2 ) + (4.59) 

iV 

then Taylor's theorem in the form of proposition 14.131 and hypothesis 12 . 71 give the upper bound 

\^ Z (<&[/ K(47 (x, k),j»(u, z)) - f{z)\) I 

\& z (c 2 N { 7 (x, k)-x}. V f (z)) I + HI VVf#\\t E Z (4l|fc|| 2 ) 

+ C 2 Ar ||V/#||#E z (||^( U ,z)-^ r (z)||lL [1 , it(/)] (7r m (z))) (4.60) 

<Lo|||/ # ||| # + |||VV/#||# f^+^supc/) +^ (/) ||V/#||#7r m (z) (4.61) 



where 



^L7r m (^)^ (/) |||/#|||# (4.62) 
L = Lo + l + ^a^+b^supc^ 2 ). (4.63) 

N 



I 



Proposition 4.16 Suppose f <G $ C) cn — > 00 as iV — * 00, Z ~ iV(0, J^) and z = (7r m (z), x, (1 

o- u n(z)Z 6 u «(z) 

i)7rr(z) + tj£b(Ug , z)) /or i € [0, 1] and fc = — 5 1 ^ — then for each s € [0, 00), 

c N c N 

di I di 

I N (z,t) = J2^ Z h(x, k)-x-k] i [d i f(z t ) + ^[i(7(^, fc) - x + ] ^0 (4.64) 

uniformly on € x [0, 1] as ./V — > 00. 

Proof. We may assume that N > No where for all N with N > No we have that ^ s + 1. Set 
(T 5 = ^(an 15 ), if z £ €. s then 

^(-M)! «S 25- 2 E z ( c 2 r ||fc|| 3 ||V/|| co + c 2 r ||fc|| 3 a||VV/|| 0O ) (4.65) 



SdlV/IU + tfHW/lloo) (* w s,„„„3 



< , 2c ; - {^(\\z\n + ^)^o (4.66) 

uniformly on (£\£ 5 ) x [0, 1]. Now consider the more tricky case when z S C 5 . Let z\ = (7r m (z), (1— 
t)7Tr(z) + tjf^(u^,z)). We drop the (z,t) dependence for notational convenience and write I N — 
7f + If + Jf where 

di I di \ 

Ii = ^c 2 N E z [ 7 (x,fc)-x-fc]49J(z t )+^[i( 7 (x,fc)-x + fc)],a 2 a j /(z)]l|| fe ||^ ; (4.67) 

= c 2 r E z ([ 7 ( 2 ;,fc)- a; -fc]-V/(zJ)a|| fe || <5 ); (4.68) 

di / di \ 

^ = [7(^A : )-x-fc]4aj(z t )-9J(zJ) + E[i( 7 (x,A : )-x + fc)] J 9^ J /(zp|| fc || < , , 

i=i V j=i / 
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(4.69) 

then Ii obeys the same bound as inequality l|4.66|l and therefore vanishes uniformly on <L S x [0, 1] in the 
limit N — > oo. If 7(2;, k) — x 2 then 7(2;, k) — x — k = — 2$ifco • n{xi)n{xi) and therefore the boundary 
condition on the class of test functions: n{x\) ■ V ' f{z\) — fdcfimtion ^.llfl implies that 

If - c 2 N E z ([-r(x,k) -x-k}- Vf(zi)l l(x ^ X2 l lM<s ) . (4.70) 

Proposition 14.81 implies that j(x, k) ^ X2 only if \n(x) ■ k\ ^ (Ao + i?o)||fc|| 2 - If 7(2;, k) ^ X2 and 
J2^=o^ n ~ ll^ll then proposition l4~7l gives an estimate for \(j(x, k) — x — k) -\7 f(z)\. Ifj(x,k) ^ x 2 and 
SJ^Lo ^ n 7^ ll^ll then we may still use proposition 14.71 and add an extra term given by proposition 
In consequence 

1(7(3, k)-x-k)- Vf(z\)\ < 1(7(3, k) - x - k) • V/(**)| + \(j(x, k)-x-k)- [Vf(z\) - V/(z*)]| 

(4.71) 

/ oo \ 3 / oc \ 2 

^ 2(2fl +^o)||VV/||oo(5^1»n) +2S (llfcll-j^On) ||V/||oo 

\n=0 / \ n=0 / 

+ 2||VV/|| co ||fc|| 2 (4.72) 
< [2{(2B + ^o)'5 + l}||VV/||,3 C + 2So||V/|| 00 ]||fc|| 2 . (4.73) 

Accordingly 

< [2(2B + A )<5||VV/|| 00 + 2i3o||V/|| 00 ]E z ( c 2 v ||A ; || 2 a|„ (;c) . fc ^ (Ao+So) || fc || 2 ) (4.74) 

< [2(2^ + A))^ W/IU + 2Bo\\Vf\\ OQ }F 2 (A Q + B ) 

CNU s+l,R(f) 

uniformly on £ 5 x [0, 1] using corollary E3 Split = if + if where 

= J2°n eZ l^ty-x-kUdifiz^-difizti+^ihti&ty-x + VWifW 
i=\ \ j=l 

x l 7 (x,fe)^ 2 l||fc||< ( 5^ ; (4-76) 

J 5 = [7(x,A;)- a ;-fc] 4 [a 4 /(z 4 )-9J(zJ)+^[i( 7 ( a :,A ; )-.T + fc)] J ^9 J /(z)] 

i=i \ j=i 

x l 7 (2;,fc)=^2 I1 -||fc||<<5 (4.77) 



and notice that 

|7f| < 4||VV/|| 00 E z ( c 2 v ||fc|| 2 a K;c) . fcK(Ao+So) || fc || 2 ) (4.78) 
4||W/|UF 2 (A + i ?0 )^ o 

LN "s+l,R(f) 
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uniformly on <£. s x [0, 1] using corollary 14. 121 Now decompose 1^ as 1^ = I§ + 1^ + 1$ + 1^ where 



f d 2 \ 
J 6 = 2 X] fiZ ||fc|2 fc ' "O^X* ~ k ' n ( a; i) n ( a; i))i n ( a; i)i^^/( z ) IL 7(2;,fc)#i :2 l||fe||<<5j 5 (4.80) 

i? = 2 p z (^pjj 5 * ' n(a;i)(fc ~ fc ' < x i)^))Axi)&djf(z)hn>f>) ; ( 4 - 81 ) 

4 W = ^E^^^fe)-^-^^/^)-^/^) 



di r 



| (7(2, fc) - a; + k) - TrjjXk - k ■ n(xi)n{xi)) 

\\K\\ 



f(x,k) = X2 



h\k\\<s 1(4.82) 



di / I di 

= J2 EZ \c 2 N h(x 1 k)-x-k] l [ dJiz 1 ) ~ fyf (z{) - x^d^fiz) 

3=1 



X H- 7 (2:,fc)=2:2 I l-||fe||<5 (4.83) 



and 



(4.84) 



We observe that l/^l is bounded by a half of the RHS of inequality (|4.79() and therefore it too tends to 
zero uniformly on <L S x [0, 1]. We bound if easily: 

\I? I < 2 ( r 1 ||VV/|| 0O E* (4||fc|| 3 ) < 8||V J, /lloC ^ (<£E*(W) + §0 - (4.85) 



A* 



-N 



uniformly on £ s x [0, 1]. Note that 



3 = 1 



f> - xi] s f [d^fisz 1 + (1 - s)4) - dtdjfiz)] ds 

3 = 1 J ° 

(4.86) 
(4.87) 



di 



< ||fc|| 5>(c^/, \\k\\+c N 2 H% f) R(f)) 
3=1 



and 



di di 

J^wididjfMl+c^H^Rif)) < Y, w ( d ^f^o \\Z\\c N 1 + b 00 c N 2 + c N 2 H^ f) R(f)) (4.88) 
3=1 3=1 

which is bounded and, by the uniform continuity of didjf, tends to zero uniformly on € s x [0, 1] at each 
Z €R dl . Therefore 



di 



\I»\ < 2E Z [clWkfY.w^d,!, a 00 \\Z\\c N 1 + b 00 c N 2 + c- N 2 H™ u) R(j)) 



(4.89) 



3=1 



di 



< AE Z \ [aU\Z\\ 2 + b 2 00 c N 2 }Y,Md l d J f,a 00 \\Z\\c N 1 +b 00 crf + c N 2 H^ f) R(f)) ] (4.90) 

3=1 
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which goes to zero uniformly on £ x [0, 1] by the dominated convergence theorem. Let 



tN 

2 10 



-2 ]T E z (<&jjfjj2* ■ n ( x )( k - k ■ n(x)n(x)) J n{x) l d l d ] f(z)^ . (4.91) 



Then proposition 14.11 implies that 

tN _ 
'10 J 9 



\I? -m < 8AoE z (^||fc|| 3 )||VV/||, 

^ — 1 f _3 T\?Z { \\ >y\\3\ i l3 — 3 



< 32A c^[^E Z (ll^ll 3 ) + Cc w 3 ]||VV/|U - (4.92) 

uniformly on C 5 x [0, 1]. We need to show that — ► uniformly, to this end set Z to be the rotation 
of Z by angle tt about the axis n(x): 

Z = 2Z-n(x)n(x)-Z, (4.93) 

and notice that for any integrable function g(Z) we have E z (g(Z)) = E z (h[g(Z) + g{Z)\). We also 
define _ 

_ <t,.n(z)Z b„N(z) 

k = v ' + -^r^, (4.94) 

c N c N 

so that k ■ n(x) — k ■ n(x) and 

_ _ 2(6 u jv (z) — b u N (z) • n(x)n(x)) 

k — k ■ n(x)n(x) — —(k — k ■ n(x)n(x)) H 5 ^ • (4.95) 

C N 

Thus 7$ = I" + 7$ + J* + l£ where 

dl ( t \ 

Iii = -2 E EZ y k ■ n ^ " M 2 ) ' n(x)n(aO)X*)«W(z)J ; (4.96) 



X I1 -fc-n(i : )<-Bo{||fc||v||fe||} 2 ) i (4-97) 



X a -B {||fe||V||fe||} 2 <fc-n(a:)<0 ) ' (4-98) 



and 



j£ = E fiZ f c ^ fc ' n W f-jj|jja " ( fc " fc ' WMlfc-nC*)^ ■ (4.99) 



j,j'=i \ 
We continue to bound: 



2& 00 ||VV/|| 00 E Z (||*||) < 2M|W/l|o ° (a^iWZW) + ^ ) ^ (4.1(H)) 

Cat \ C N J 

uniformly on € s x [0, 1] as N — > oo. Furthermore if fe • n(x) < — _Bo||fc|| 2 then for all r £ [0, 1] we have 
that 

dl .1 

w(a; + rfc) = + rk ■ Vw + E \ C 1 _ s)k i k j d i d : jUj{x + srk) ds (4.101) 

. , JO 

«j=i 

«S w(a;) -r{||Vw(x)||B - ^HVVwIloo} ||fc|| 2 (4.102) 
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where we have recalled that B$ — |||VVw||oo and ||Vo;(x)|| ^ 1 as z € C 5 . It follows that $o = ||&|| 
and therefore that $1 = 0. This and the corresponding remark for $i imply that = 0. Turning our 

N 
13 



attention to I?i we derive the bound 



< llVV/IUE^^pfli^i^iifeiiO + llVV/IUE^^IIfell 2 !,^^^,^,) (4.103) 

+ 2 " W {"^. (4.104) 



c -Af (7 s+ l i _R(/) 



-AT 



This tends to zero uniformly on € s x [0, 1] by using proposition 14. Ill 
Next we consider the expression 



0? K 



If #i , #i ^ we deduce that 



—2 



#f #1 



\\k\\ 2 \\k\\ 



k ■ n(x) 



\W \\k\\ 



k ■ n(x). 



(4.105) 



\\Vu(x)\\ 



\\k\\ 



di 



h y 



jrrrhkj / (1 — s)didjLu(x + S<&ok) di 



$ 2 __ /■! 

_° fcjfcj / (1 — s)didjLu{x + sfiok) ds 



< ||VVa;||oc(||fc|r + ||fc|| 2 ) < 2(A + B )(||fc|| + \\k\\y, 



(4.106) 
(4.107) 



whereas if t?i = and #i ^ so that u)(x + i^o^o) = and ??i = 0, we may set x\ = x + t^o^o and 
if = k — t?ofco to derive 



■jj-Tu/s ' "■( a; ) — & ' n(x) — -rrrrrk ■ n(x) 
\\k\\ \\k\\ 

— K ■ n(xi) + K ■ (n(x) — n(xi)) 

< A \\K\\\\k\\+B \\K\\ 2 ^ (A Q + B Q )(\\k\\ + \\k\\) 2 

where we have used proposition 14. II and the result that if K ■ n{x\) > Bo\\K\\ 2 then 

uj(x + k) = uj(xi + K) > oj(xx) + K ■ Vu){x\) - B Q \\K\\ 2 > 0. 

This allows us to conclude that uj(x + k) ^ f2. Accordingly if k ■ n(x) > then 



0\ 



WW Hfcll 2 / 
In consequence we are led to a bound for Iy± 



k ■ n{x) 



^2{A^B G )(\\k\\ + \\k\\f 



(4.108) 

(4.109) 
(4.110) 



(4.111) 



(4.112) 



N . ^ 8(A)+5 )||W/|| 



00 E z 



cn 



<Too\\Z\\ + — 
CN 



uniformly on € s x [0, 1]. 



(4.113) 
I 
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Theorem 4.17 J//6j c U {/ } wifA extension f* e C 2 (C # ), itf ->• s and cjv — * oo as N — » oo t/ien 
lim E z (<&[/ (7r m (z),7 fc) ,3&(u?,z)) - f{z)\) 

= ±a s (z)Af(z) + b s {z) ■ V/(z) + H s (z) • V/(z). (4.114) 
<j u n(z)Z b u N(z) 



N- 



uniformly on € where x = tt-^(z), Z ~ N(0, Idi) an d k 



c N 



"JV 



Proof. The case / = /o is trivial, so consider / € 5c- By remark 14.141 we must prove that 
E z (c^fG N (u^ , z, Z)) — ► uniformly on £ as N — ► oo. In terms of the modulus of continuity of 

^i/ w e compute 



«)-*!"(*)]• j n E z (^f*(n m (z),(l-t)x + t 1 (x,k),(l~t)n r (z)+tj^(u,z))-Vf(z)) dt 

*oo\\Z\\ , ^ , i^ (/) i?(/)\\ 



-JV 



uniformly on £ by noting that 



V # , + ¥ + < 2|||V/#||| 



(4.115) 



(4.116) 



and using the dominated convergence theorem together with the fact that dif& is uniformly continuous. 
Similarly, 



e 2 N E z (k) • / (V/(7r m (z), x, (1 - t)Trr(z) + z)) - V/(z)) dt 



&»(*)• / (Vf(TT m (z),x,Mz){l-t)+3^{u,z)t)- V/(z)) dt 
H R(f) R (f)\ „ 



i=i V 



(4.117) 
(4.118) 



uniformly on £ as jV — > oo, 

'0 «/o 



E z (( 7 (a;,/c)-a;) J (7(a;,/c)~a;) J 



x (didjf*(ir m (z), (1 - t)z + t 7 (ar, fc), (1 - s)^(^) + sjg,(u, z)) - c^/(z))) dt ds 



di 



# 



>IIZII , fe^ , #H (/ )i?(/r 



cat c w 



"JV 



and 



2 \ - bjbjdjdjfjz) 
c N 



E 



2c' 



2r 2 
ZC JV 



(4.119) 



(4.120) 



uniformly on £ as jV — > oo. Using proposition 14 . 1 0l it remains only to prove that 



c 2 N / E z ([ 7 (a;,fc)-s-A : ]-V/(^ m (z),2;,^rW(l-t)+JiDK^)) dt 
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+ \(? N ^ Z ((l(x,k)-x-kU 1 {x,k)-x + k) 1 d l d J f(z)) ->0 (4.121) 
uniformly on £ as TV — > oo which is a direct consequence of proposition 14.161 



5. Convergence of the Sequence of Stochastic Processes 

In this section we use a relative compactness argument to deduce the existence of a sequence of 
measures that converge weakly to a weak solution to equation l|l.l|l in the case of the Linear Growth 
Condition or a vague solution to equation (|l.lfl if the more general growth conditions are satisfied. 
The strategy is to consider spaces of discrete measures with bounded total mass and investigate the 
two component process (P N ,u N ) where ti" is an approximation to the time. By using u N instead 
of the true time we may utilize a stochastic generator that has no explicit time dependence and it is 
moderately easy to write down a set of rules that solve the associated martingale problem. We proceed 
to show that the sequence of stochastic processes is tight and Prohorov's theorem reveals the existence 
of an accumulation point. We verify that this accumulation point is a solution to equation l|l.l|l in the 
relevant sense. Once this is done we prove a result that allows us to dispense with the fictitious time 
variable. 

It is advantageous to introduce an extra process associated with the fictitious time variable, the 
purpose of which is to bound the total rate of all processes away from zero. This clock process does not 
alter the current state of the measure describing the stochastic particle ensemble, it merely advances 
the fictitious time variable. In this way we can be sure that if all the other rates momentarily become 
zero the computed holding time in the current state remains finite. 

Throughout this section we assume that tun and cjv are sequences of positive real numbers that 
tend to infinity as N — > oo. 



Definition 5.1 Define spaces of measures by 
1 C = {P£ M B {£) ■ P(vr m ) < C}, and 



m N = |p = J25 Zi : Zi e £, P(7r m ) m N j 

(5.1) 



where S z is the Dirac measure concentrated at z. It is convenient to set E ~ 971 x [0,oo) 



Definition 5.2 We put the metric D(-,-) on E = Mb(£) x [0, oo) and its subspaces in terms of 
equation H2.24JI by defining 



D((P,u),(Q,v)) = 

"'weak 

(P, Q) + \u - v\ 



(5.2) 



and let D be the Skorohod metric on the space of cadldg functions De[0, oo) derived from D where if p 
is a metric on a space X then the Skorohod metric p on Dx[0,oo) is defined by (cf. equation (3.5.2) in 
1) 



p(x, y) = inf 
AeA 



7(A) V / e sup p(x t at, y \(t) at) dT 
o t^o 



7(A) = ess sup 

0^t<s 



log 



A(s) - A(t) 



s - t 



(5.3) 



where 



A = {A € C([0,oo), [0,oo)) : A is a strictly increasing Lipschitz function with 7(A) < 00}. 



(5.4) 
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Remark 5.3 We adopt the notational conventions for the mass-weighted source moments: 
A t = I t (t) and A<«> = sup{J t (7r&>) : < t < oo}. 



(5.5) 



Definition 5.4 Let v™ : M B {£) X 5B(C) n -► R oe denned byu^{P,Bt x •• • xB„) = N- n u n {NP,B 1 



x B„) where 



u n {P, Bx x ■ ■ ■ x B n ) = CO) H n . . . n 



(5.6) 



disjoint cycles 

(il,...,tfc) of CT 



Pere S n is the symmetric group on { 1 , . . . , n} and £ is the signature homomorphism. 
Remark 5.5 By proposition lC.il 



v i ( ]v 51 *«* ' ^ X " ' X ^ J = 



fc=i 



?y— t if there exists an injection a : {1, . . . , i} — > {1, . . . , n} 
x such that z r = w a ( r ) for r = 1, . . . , i; 

otherwise. 

(5.7) 



Definition 5.6 Define the sequence of generators A N : M(E N ) — ► M{E N ) to be 
A N (F)(P,u) = Nr{F(J^ K (P,u))-F(P,u)} 



+ N / {F(jf ({P,u),z))-F(P,u)}I u (dz) 



T J 



+ Ncff / E < F Jq (P, u), z 



cat c|, 



F(P,u) mdz 



»=1 3=1 

x {F(J%((P, u), zi, . . . , zi, wx, ■ ■ ■ , Wj)) - P(P, u)}^(P, dz u . . . , dz,) (5.8) 



w/iere Z - N(0,I dl ), and the jumps: J£ LK : F — > F, : F x £ — > P, J^E^xCx 
and PX : £ w x (P x <V — > P /or i = and j = 1, . . . , J are given by 



J?{{P,u),z) 



J»((P,u),z,k) = 



J i!j((P; u),Zl, ...,Zi,Wl,..., Wj) 



P,u 



1 



P N (Pu) 



1 



P N (Pu) 

P(7T m ) ^ 

otherwise; 



P 



p N (P,u) 



(5.9) 
(5.10) 

(5.11) 

; (5.12) 



p+Uf^-i^X 

n—1 / 



l 



p N (P,u) J 
fori 1 /../ 1 J/: (5.13) 
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and 

p N (P,u) = N |r + A u + 4P(£)+^^^ K Uti , j (z 1 ,...,z i ,€P)^(P,dz 1 ,...,dz i ) | . (5.14) 

Notice that the image of each of the jumps is a member of E whenever P({z}) ^ in equation (I5.12|) 
and (P, {zi} x • • ■ x {zi}) ^ with J2m=i Ir o(%) = ELi 7r m(zn) w equation <|5.13[1 . TTiis provides 
the motivation for definition \b~.l\ 

Remark 5.7 The generators define a standard Markov jump process: 

A N (F)(P,u)=p N (P,u) f {F(X)-F(P,u)}fx N ((P,u),dX) (5.15) 

J E N 

with transition function p N : E N x %$(E N ) — > [0, 1] defined by 

rN 

fx ((P,«),-) = pJ v (P)U) ^ff LK (P, u ) 
'at /" 

+ ^(p^)i c 7 ™ (d ^((^) 



p^(P, W ) 



P( 



TV 1 J 



>(dz)E z (d^ ((Pjtt) 

pW( P 7 u) E E i!7T W*. • • • , dtBx, . . . , dlBj) 



x^f (P dzi, . . . ,d^)5jjv .((F,„), zl ,... A , Wl ,... >t „ 3 .). (5.16) 



Note that 



p N (P,u) ^ N{r + + c 2 N m N +lJK oc m I N + K™ N ). (5.17) 

In consequence the generators are bounded for each N and according to section 4.2 23^ there exists a 
sequence of i^-valued random variables (P t Ar ,u^) such that 

Mf (P ) = P(Pf ) - P(P JV ,<) - /V(P)(PjVf )ds (5.18) 

is a martingale for each F S M(E N ) for any random initial condition (P^, ) € P w . 

Hypothesis 5.8 Let Uq = /or all N, and suppose that there exists a sequence P^ € 9Jl N with 
P^ — >■ Pq as N ^ oo such that P^ is deterministic with P^ (Tr m ) bounded above. For notational 
convenience we define 

S = supP JV (^ m ). (5.19) 

AT 

Remark 5.9 By considering the bounded measurable functions tt^(P,u) =uAL and $/(P, u) = P(/) 
where f £ $ we find that 

M t "(7tf ) =v?AL-J* p N (P s N , u?) | (V + -pripryrj ) A i - «f A l| d s (5.20) 

and 

Mf^/) = P t N (f)-P N (f)- [ [ f(z)l P N iVm)+nAz)/N<m „I u? (dz)d S 

JO J£ 
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-c% [ f E z {/(^ m (z), 7 ( 7 r n (z), CTu (z)Z Cw 1 + 64z) Cw 2 ),j£(uf,z))-/(z)}Pi v (dz)d S 

JO J€ 

~H~T1 i i K u N^ j (z 1 ,...,Z l ,dw 1 ,...,dw j )\ f(Wm)-J2f{Zn) \ 

i= i i=i J o ^C*+j [ m=1 n=1 J 

x i/f r (Pj v ,d«i,...,d8 i )da(5.21) 

where Z <~ Af(0, J^), are martingales for each AT. 

Proposition 5.10 Let (J r £ r ) t>Q &e a filtration adapted to (P N ,u N ). Suppose that ir u (P,u) = u and 
7ri 2 ^ (P, u) = u 2 then 



M t N (n u ) = u?-t and M t N (^) 



(5.22) 



are -martingales for each N. 



Proof. We have 



E(t^AL) = t+ f t E(A N (^)(P T N ,u^)-i)dr 
Jo 



I 

Jo 



t+ I Elp N (P T N ,u?)i [< + 



.N 



1 



P N (P T N ,u?) 



A L 



1 



(5.23) 



dr < t. (5.24) 



p N (P T N ,u?) 

Hence by the monotone convergence theorem is of class L 1 and E(u^) ^ i. Furthermore if t ^ s, 

1 



E|E (uf A L - t | J^f) - (uf A L - s 



E 



Elp^,^ 



1 



;|(t-«)p ^ 



^(p*,<) 



•ff ) dr 



- ! > («, v >/.- — 

< (t-s)E(E(u?\^j) < i(t-*) 



L-N~ l r- 







(5.25) 
(5.26) 
(5.27) 



as L — > oo, the conditional form of the monotone convergence theorem and Fatou's lemma then imply 
f - t is an ^ 



that — t is an .F^-martingale. Similarly 



E ( (uf) 2 AX 2 ) = ^ + ^E^ ( P-n- ) |(< + -^^ 



AT) 

(uf) 2 A L 2 



A L 



2u 



iV 



dr (5.28) 
(5.29) 



Thus the monotone convergence theorem gives that E((uf^) 2 ) < oo. Moreover 



E 



E[(u?) 2 AL 



p N (P T N ,u?){(u? + 



P N (P T N ,U?) 



(u?) 2 }dr\^ \-KYal 



-N 



N\2 
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< (t - S )E {E ({2< + N-'r- 1 }!^^^ \ F») } (5.30) 
= (t - s)E ({2uf + iV-V- 1 }!^^^^) - (5.31) 

as L — > oo by the dominated convergence theorem. Applying the conditional form of the monotone 
convergence theorem and Fatou's lemma on the LHS we deduce that 

E (V) 2 - j\ N {P?,u») | (V + pjV(p ^ <} ) 2 - «) 2 j dr | ^ = (uf ) 2 a.s. (5.32) 

and thus M i JV (7ri 2 ' ) ) is an JF^-martingale. I 

Proposition 5.11 For aZZ TV and i ^ 0, 

E((P 4 JV ( 7 r m ))' i ) < Qi(t) » = 0,...,J, (5.33) 
where Qi is an increasing polynomial on [0, oo) of degree less than or egual to i. 

Proof. We prove the result by induction. If i = 0, then we can take Qo(t) = 1. For < i ^ I, 
we assume the result holds for all r < i, set F((P,u)) — (P(ir m )Y < m % N in equation (|5.18|l and take 
expectations to deduce that 



E((P/> m )y) = (i^(7r m ))* + JV 



o 



A 

x ^-P^{n m )+-K Ta (z)/N^m N 4f (dz) ds (5.34) 

< g + E f ;„ rVr i / QrW^^Q^t), (5.35) 

r=0 * )■ ■ JO 

where we have made use of hypothesis 12.51 I 

Proposition 5.12 For all T > 0, / G $ c U {/o} i/iere exist an increasing polynomial Q on [0, oo) of 
degree no more than 2+1 such that 

(a) E^upJMf^)! < (^Jp-) 2 KR(f) H R(f)\\\f*\W*-> ( 5 - 36 ) 

(b) E sup |M t ($ OTin )| < 2(^ r ) (5.37) 



(c) for all e > i/iere exists Rt{^) > suc/i £/ia£ 



sup |uf - t\ > R T (e) ) ^ e. (5.38) 



(d) In particular 



sup |uf -i| ^ AT* ) <2(T+l)*r" = iV-i -► as A ^ oo. (5.39) 
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Proof. The previsible increasing processes associated with the martingales equations I|5.20I5.21|) are 

given by 



p-(Pj\<) / {F(X)-F(P s N ,u?)} 2 ^((P»,u»),dX)ds 

E N 



(5.40) 



= rN f{F{J^{P?,^)) F(P?,u? )} 2 ds 
Jo 

l z {F(J»((P a N ,u?), z, a u » (z)Zcr N 1 + b u ,(z)c- N 2 )) - F(P S N , u^)} 2 P s N (dz) ds 



Nc 



N 



J£ 

{F(J?((P s N ,u?), z)) - F(P S N , u? )} 2 4» (ds) d s 



to J<t 
i J 



K u N ti j(zi, Zi,dwi, dwj) 



x{F(J"((P s n ,u^),z 1 ,...,z u w 1 



, . . . , UJj 



)) - F(P? , )} V (P S N , dz u . . . , d*) ds 



(5.41) 



leading to 



(M JV ($ / )) t < 



f{z) 2 I u ,{dz)ds 



N Jo 

i J 



E z {/(7r m (z), 7(7i^-(z)> ( z ) Zc n + ( z ) c at 2 ) > iro(«f > *)) 



/(z)} 2 Pf(d^)d S 



,-(2 1 ,...,;z j ,dt0i,...,dw ; -) 



E/( w ™)-E/( z ») ^f(ff,d^,...,d Zi )d S . (5.42) 



n=l 



Write fc = cr u jv (z)^^ 1 + b u N(z)c N 2 and a; = 7Tq(z). The inequality 

cf v E z ([/(7r m (z), 7 (z, *)) ~ /(*)]') 

^2 C 2 v E z ([/( 7 r m (z), 7 ( 2 :,fc),jS(uf,z))-/(^ m (z), 2 :,jS(uf,z))] 2 ) 

+ 2 C 2 r [/( 7 r m (z), a; ,^( U f,z))~/(z)] 2 

^2c 2 v E z (||fc|| 2 ||V/#||# 2 ) +4||/|| 00 H^ /) ||V/|| 00 7r m (z) 



s$4 



then shows that 



t&dt + foLsupc^ 2 ) ||V/ # ||# 2 + ||/|U^ (/) ||V/|| 

S(t) 



(M ($/))t ^ 



4iV 



(5.43) 
(5.44) 

(5.45) 



(5.46) 



where 
S(t) 



AA^t 



16^1 + 166^ sup c A 2 + 16 



N 



Jo 

+ 4 j* Pj> m ) ds + 4ZJ(1 + jfK x f % (P s N (n m f ds} K% f) H% f) \\\f#\\\* 2 . (5.47) 
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It follows by proposition 15.111 that if 



Q(f) = 4 



4o*,di+4i& supc 2 + 5 



/ Qi{s)ds + AA {a) t + 41J{l + jfK^ f Qj(s)ds (5.48) 
Jo Jo 



then E(Af Ar ($ / )) T < jQ^N' 1 K% {f) H™ u) \\\f*\\\* 2 . Doob's ^-Inequality 



E sup \M?(F)\) ^4E(M N (F)} T 

O^t^T 



(5.49) 



gives part (a). Our definitions R(fo) = and — = 1 are essential to check the n = case. 
Furthermore inequality l|5.42|l and the mass preserving property of the interaction kernels given by 
hypothesis EHI give (M N (®„ m )) t ^ tA^/N, and thus inequality (fB"^ leads to part (b). 

By proposition ^. 1UI (-k u ) and M^(tt^) are martingales. In consequence the increasing process 
associated with M^{-k u ) is given by equation l|5.41|l with F — tt u , thus 



rl\ r 



and hence 



E sup | u f-,sK2( — 



rN 



(5.50) 



(5.51) 



Inequality (|5.38() follows from inequality (|5.49|) applied to 



t and the Markov inequality, with 



i?r(e) = 2T?r 2 £ 1 establishing part (c). If we set e = 2(T + 1) 2 r ?N a then part (d) follows. 



Proposition 5.13 There exist an increasing polynomial Q on [0, 00) of degree no more than X such 
that (a) 

for all f € 3c U {/o} t cmd ('fej 



i/($ Im )(^,«f)ia (1) 

/or a// f ^ 0. 

Proof, (a) Setting x = 7Cq(z), k = <j u n {z)Zc n 1 + 6 u «(z)c^ 2 and Z ~ N^O,!^) we have 

(^(^(Pf^f)! < /|E z (/(^ m (z), 7 (x,A ; ),^( U f,z))-/(z))|P t w (dz) 

.ij^i) • • • >2i>dwi, ■ • • ,dwj) 



(5.53) 



1 J 



f( w m) - 



rn — 1 



v?(P t N ) dz 1 ,...,dz i ) (5.54) 



< A( )||/|| oo + Lfl^ /) |||/#|||#J^(7r m ) 

TO'OO dJV/_ Ml j-II 1 T S7 (' 



+ k™ u) p? (7r m )||/|U +x^r(T + j-j^jf (tt^h/i 



(5.55) 



where we have used proposition 14.151 with L given by equation l|4.63[) . Accordingly equation (|5.52|) 
follows with 

Q(t) = AW + (L + l)Q 1 (t)+lJ(T + J)K 00 Qi(t) (5.56) 
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where we have made use of proposition 15.111 The n = case is dealt with by using the definitions 
H$° = K%° = 1. 

(b) Inequality (|5.53|) follows directly from inequality (|5.54|) together with the mass conserving property 
of the interaction kernel (hypothesis I2.8|l . I 

Proposition 5.14 For all T > and e > there exists such that for all N ^ 1 we have that 

supP t N (7T B1 )> L T (e)) <e. (5.57) 

Proof. If t ^ T we have from the martingale definition 

Pf(7r m )^P JV (7r m )+ [ T \A N (^J(P t N ,u?)\dt + \M t N (<^J| < 5 + TA«+sup|M^($OI- (5.58) 
Take the supremum and expectation and apply proposition 15. 12T b) to find 

supP t N (TT m )) ^K T (5.59) 
where 

K T = S + TA« +2(TA^)i (5.60) 
and thus, if Lt(e) = Kt/c then Markov's inequality gives equation (|5.57|l . 1 

Corollary 5.15 For all T > and < n < 1, 

infP((P t Ar , U f)eOT iT(l;/2) x [0,r + i? T (7?/2)], 0<t<T)^l-7?. (5.61) 

Proof. This is immediate from propositions 15.121 and 15. lH 1 
Proposition 5.16 For all C > 0, the space 9Jlc is a compact subset of Mb(£)- 

Proof. First we show that 9Jlc is closed. Suppose p n E Tic with p n — p, then for each L 6 N set 
7r m( z ) = L h 7r m (z), we have that C > p n {^^) — > p( 7r m)- Thus p{~K^) ^ C and by taking L — * oo, the 
monotone convergence theorem implies that p(7r m ) ^ C, i.e., p G 9Jtc since p(£) p(7r ro ) $5 C and so 
pe Mb(£). 

Next we show sequential compactness. Suppose p n G 5D?c> write p„ = A„g„ with q n (ft) = 1 and 
A n = Pn(£)- Then A„ 6 [0, C] and therefore there exists a convergent subsequence A rlfc — ► A. If A = 0, 
then for any bounded continuous function / on £ we have Pn k (f) ^ A n J|/||oo —> as k — > oo, i.e., 
Pn fc — 51 0. If A 7^ then we may assume, by passing to a further subsequence if necessary, that X Uk > ^A 
for all k and we find that for any rj > we have 

q nk (tt-^I, • [2C/A771}) = 1 - q rik (v' 1 {[2C / \ V ] + 1, ...})> 1 - ^Cg„,(7r m ) > 1 - i, (5.62) 

where |~x] = min{n G Z : ra 5^ x}. Note that since 7r m is continuous 7T~ {1, . . . , |~2C/A?7]} is closed 
and bounded, and thus compact. Accordingly the g nfc form a tight sequence of probability measures. 
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Prohorov's theorem (theorem 3.2.2, [231 ) implies that there exists a further subsequence 9n fc(r) which 
converges weakly to a probability measure q. Hence 



dwcak(p™ fcM , Xq) = ^2 ^(|An fc(r3 «n fc(r) (/0 - A 1) -> as r -> oo (5.63) 

with limiting measure Xq € 97lc by closure. It follows that 9Jlc is sequentially compact. | 



Corollary 5.17 The sequence (Pf,u^) obeys the compact support condition (remark 3.7.3, \23^). 

Definition 5.18 The modulus of continuity of an element x € De[0, oo) is defined to be 

w (%, 5, T) = inf max sup D(x s , xt) (5.64) 

'« i s,t6[t«-i,*0 

where — to ^ ti ^ ■ ■ ■ ^ t n -\ ^T ^t n and min{/j; — > <5. 

Proposition 5.19 Por any T > and 77 > 0, 

]hasapF(w'((P N ,u N ),5,T) > ?y) 77 (5.65) 

JV 

/or sufficiently small S. 

Proof. We investigate the partition given by tj = 2j5, then for x £ De[0, 00), 

w'ix, 5, T) ^ 2max sup D(x2j8,x s ). (5.66) 

J 2j5<s<2(j+1)5 

It follows that 

F(w'((P N ,u N ),6,T)>ri) < P(max sup D((P 2 *, (Pf, uf )) > ±77 ] (5.67) 

\ J 2j«<s<2(j+1)« J 

^ -Ejmax sup D((P 2 N js ,u» 5 ), (P s N ,uf )) j . (5.68) 

1 y J 2j<5^s<2(j + l)5 / 

If 2j'5 T and 2jS < s < 2(j + 1)<5 and we write $/„ (P, zt) = P(f n ) we have 

< J2^(\ P s N (fn)- P 2 N j s(fn)\^^) + K ~<S\ (5-69) 



4 

n=0 



oo 



n=0 



Mf(d> / „)-M^($ / J+ / ^(W,<)dT 

2j<5 



A 1 



+ k"-<- 5 | (5.70) 



n=0 



1 ( 2 sup |M T " (* /n )| + f |^($ / J(P T ^,<)|d, 



A 1 



Os£t^T+2(5 

Thus propositions 15.121 and 15.131 give 



I2j5 

+ 2 sup \u?-t\ + 26. (5.71) 



E max sup D((P 2 ^ , u&), (Pj v , <)) 

V J 2j«<s<2(j+l)« . 
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oo 1 

V- 

Z-^i An 



71=0 



2E sup |Mf ($ / J| + 25Q(T + 25)^ (/n) if~ (/n)llu „ 



f#\\\# 



A 1 



■ 2E sup \u? -t\ + 25 



(5.72) 



E- 

L^i An 



71=0 



2K^ fn) H^ f J\f*\\\* 



Q(T + 2(5) 
TV 

<5Q(T + 25) 



1 ' 4 ' 2 +2S (573) 



where we have made use of inequalities (15.36(1 and ((5.50(1 and checked the n = case separately. 
Consequently 



limsupP(u/ \(P N ,u N ),6,T) 



4<S(l + g(T + 2)i;~ 2-* n ) 



N 



whenever 6 ^ W[l + \/2Q(T + 2)/(\/2 — l)] -1 A 1. 



(5.74) 
I 



Corollary 5.20 TTie Zaios o/ Z/ie random variables (P N ,u N ) are relatively compact in the space of 
probability measures defined on De[0,oo) equipped with the weak topology. 

Proof. Note that since £ is polish so too is E. The result follows by applying proposition 3.7.4 of |23j . 
together with corollary 15 . 1 71 and proDOsition l5.19l I 

Corollary 5.21 There exists a subsequence (P n ,u n ) converging in distribution (which we denote by 
=$■) to (P,u) £ Ce[0, oo) almost surely. 



Proof. Corollarv l5 .201 implies the existence of a subsequence converging to an element of De[0, oo). By 
theorem 3.10.2 [23 to prove the corollary it will suffice to check that the distance between neighbouring 
states vanishes uniformly as N — > oo. We have for any (p, v) £ E , z, Z\, ... , Zx, W\, . . . , Wj £ €. and 
k £ K dl that by using \\f n \\oo < |||/n||| < \\\ft \\* < 2" for all n £ N that 

D((p,v),J^ K (p 7 v)) 



1 1 

p N (p, v) ^ WV : 

fn(z) 



N 



A 1 + 



1 2 



p N {p,v) N rN' 



(5.75) 
(5.76) 



D(^v),J^({p,v),z,k)) 



E4- 

n=l 
2 1 

N + rN' 

D ((P> V )> J i!j((P> v )> z i> • ■ • wi, . . .,Wj)) 



/„(7T m (z), 7(^(2), fc), jfciv, Z)) - f n {z) 



A 1 + 



E 4 " 

71 = 



Em=l /»W " E?=l fn( z l) 



N 



N J ' p N (p,v) 

(5.77) 

(5.78) 
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Proposition 5.22 For all t ^ 0, 



sup P s (7r m ) < oo a.s. (5.80) 



Proof. Observe that Pj v "(7r m AL) =>■ P s (ir m /\L) by weak convergence, moreover, since by corollarv l5.21l 
the limit is almost surely continuous in time the convergence is almost surely uniform on compact sets 
(theorem 3.10.1 in (2H|)- Thus sup P s JV ™(7T m A £) =>■ sup P s (7r m A L). Now Fatou's lemma implies 

that 

E ( sup P s (n m A L) ] < liminfE ( sup P* n {-K m AX) J < S + tA (1) + 2(tA (2) )i (5.81) 
Taking L — > oo and using the monotone convergence theorem yields 

E ( sup P s (7r m ) ) 3 + iA (1) + 2(iA (2) )i (5.82) 

Accordingly P I sup P s (7r m ) < oo I = 1. I 

Definition 5.23 Let Vt — {(p, v) G De[0,oo) : sup p s (7r m ) < oo}, then vrovosition \5.2 e A imvlies that 
P(P G Vt) = 1 /or a// 1 > 0. 

Remark 5.24 Since £ is separable, so too is De[0,oo) by theorem 3.5.6 of j2H|. Now the Skorohod 
representation theorem (theorem 3.1.8 Ibid.) implies that there exists a common probability space 

(n,F,P) on which (P Nn ,u Nn ) are defined and (P N " , /" ) (vj) (P,u)(zu) as n -> oo for all ro G f2. 

Proposition 5.25 ///or a// 1 ^ 0, (p™, u n ) G Vt n D B w„ [0, oo) and (p, u) G C E [0, oo) n V t w'i/i p" ^ Ps 
and sup \v™ — s\ — > as n — > oo then for every f G 5c, (VJ 

(d^)ds^ / J s (/)ds (5.83) 

J£ JO 

as n — > oo and (b ) 

I \ K v n ^ 3 (z 1 ,...,z l ,dwi,...,dw j )< V f(w m ) - V f(z n ) \u? n (p™,dzi,...,dzi)ds 



n=l 

/ / i^ s ,i.j(zi, . . . ,Zi,dwi, . . . ,dwj) 

JO Jp+i 

X { ^( Wm ) ~ X! ^^™) r ^( dZl ) • • -Ps( d2 i) ds ( 5 - 84 ) 



v m—1 



as n — > oo . 



Proof, (a) For all compactly supported functions /, 

/(■2)H-pj(7r m )+7r nl (z)/iV n >mjv„ = (5.85) 
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for large enough n, as by Fatou's lemma we have limsupp™(7r m AL) ^ Ps(^m AL) ^ Ps(^m) an d thus 

n 

the monotone convergence theorem implies that limsupp™(7r m ) ^ p s (ir m ). Equation (|5.85(l follows by 

n 

the condition mjv n — > oo. Continuity of the map t h- > I t ( hypothesis 12. 5fl implies that Iv n (f) —> Is(f) 



and thus 



(dz) -/.(/) 



(dz) 



as n — > oo. We know that 



(dz) 



convergence theorem implies that 



+ |/ o »(/)-J,(/)|-,0 (5.86) 
^ Aw ll/H oo and hence the dominated 



f(z)t p n^„ tn - )+nta ^ z y Nn ^ mNn I v n(dz) - I s (f) 



ds —>■ as 



(5.87) 



(b) We make separate arguments depending on whether we are considering the self-interaction kernel 
corresponding to = (1, 1) or not. If ^ (1, 1) we have by hypothesis 12.81 that the function 



)C(s, zi,...,Zi)= / K Sti j(zi, . .. , Zi,dwi,. . . ,dwj) < f{w m ) - ^ f(z n ) \ 
J & lm=l n=l J 



(5.88) 



is bounded in modulus by (i + i)i!i!if o||/||oo'7!'m(£i) " * """m(^j)- We assume that n > no where 
sup sup u" < t + 1. Moreover since £ is separable, A / t_B(£ r ) = ®i=i -Mb(£) an d P™ — v Ps implies 

n>no O^s^t 

that (p")® r p s ® r . Then if / G # c and r < i and a : {1, . . . ,i} {1, . . . , r} we find 



1 



< ^±f^||/||oc#oci*(/r r f sup P J(7T m ) 



as n — * oo. Recalling that 



v? - (p» Bix-xB.,) = jv £ II ^np?(Bix n • • • n 

<rGS„ 



(5.89) 



(5.90) 



disjoint cycles 

Of (7 



it follows that 



/C«,zi, . . . , Zi)^"(p",dzi, . . . ,dzi) - \ K,(v", zi, Zi)p"(dzi) • --p"(d^ 
Furthermore we have 

/C(w",zi . . . , z^p" (dzi) • • -K( d2; i) - / ^( s > z 1j ■ ■ -^ijPsidzi) ■ ■ -p s {dzi) 

< / «;(x:,K-«|)pj(d»i)---p?(dai) 



/C(s,zi, . . .,Zi)p^(dzi) ■ ■ •Pg(cbi) 



0. (5.91) 
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/C(s, zi,...,Zi) p s (dzi) ■ ■■p s (dz i ) 



(5.92) 



The first term on the LHS vanishes as n — > oo by the uniform continuity of K, (which follows by the 
uniform continuity given by hypothesis 12.8(1 . equation 1(2.26(1 . combined with the fact that p"(£) — > 
p s (£) < co, the second term vanishes by the aforementioned weak convergence of (p™)® 1 . Accordingly 
for each s, 



/C(s,zi, . . . ,Zi)p s (dzi) ■ ■ -p s (dzi) 



(5.93) 



as n — > oo. The LHS is cadlag and therefore integrable in time, for large enough n the LHS has modulus 
bounded above by 



Zi)p"(dzi) • • •pj(dzi) 



< H/ll oo-^oo (1+ sup p 5 (7r m ) 1 ) <oo (5.94) 



then the dominated convergence theorem yields equation 1(5.84(1 . 

The proof for the self-interaction kernel is similar, though in this case we remark that K, is supported 
on [0, oo) x supp/ and therefore the uniform continuity of K, comes from the mass preserving property 
of the self-interaction kernel. I 

Definition 5.26 Let (P, u) be as in corollary \5.21\ then for all f € 5" c (or for all f E $ if the linear 
growth condition is satisfied), define 

M t (f) = Pt(f) - Po(f) - J W) ds ~ J Ps (K A / + b * ■ W) ds - J P s (H s ■ V/) ds 



i J 



K s>i j(zi, . . .,Zi,dw\, . . .,dwj) 



Yl f(w m )-J2f(zn) )P s (dz 1 )...P s (dz i )ds. (5.95) 



Proposition 5.27 For all f 6 # c and all t ^ 0, 

Mf n ($/) M t (/) as n -> oo . 



(5.96) 



Proof. Since (P N ™ ,u Nn ) converges in distribution to (P, u) with P((P, u) S (7e[0, oo)) = 1 we may 
prove this result by considering those sequences of sample paths with limit in Ce[0, oo). 

By corollarv l5 . 211 and the convergence in distribution of ti"" => u, an almost surely continuous limit 
in time, theorem 3.10.1 in j2H) implies that the convergence is almost surely uniform on the compact 
set [0,t], i.e., for almost all of the sample paths 



where for convenience we have written v" = u ™ (vj) 



sup \v s - s\ 

n = ,JV„ / 



as n — > oo 



(5.97) 



Let p n s = P™ "(to), with p n s ^p s . If k 



a v »(z)Z b v n(z) 



CN„ 



c 2 Nn I E z (f(n m (z)M«n(z),k)J%>(v?,z)) - f(z))p?(dz) 



with Z ~ N(0, I dl ) then 
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$a„(z)Af(z) + b s {z) ■ V/(z) + H s (z) ■ V/(z))p s (dz) (5.98) 

c%E z (f(n m (z),^a(^k),j^K,z)) - /(*)) 

- (±a fl (z)A/(z) + b s (z) ■ V/(z) + H s (z) ■ V/(z)) supp?(C) 

oo n 

(±a s (z)Af(z) + b s (z) ■ V/(z) + H s (z) ■ V/(z))^(dz) 

- / (±a s (z)Af(z) + b s (z) ■ V/(z) + H s (z) ■ Vf(z)) Ps (dz) . (5.99) 

The first term on the RHS of inequality (|5.99l) converges to zero by theorem 14 . 1 71 and the weak conver- 
gence implying p™(£) — ► p s (£)- The second term tends to zero by the weak convergence p™ — p s since 
the integrand is a bounded continuous function on C The cadlag property implies that 



/ E^(/(^ m (z), 7 (^(z),fc), Jl ^«,z))-/(z))p^(dz) 



(5.100) 



is integrable in time, and proposition 14.151 implies it is bounded in modulus by Lif^^|||/#|||#(1 
Ps(""m)) for large enough n, and thus the dominated convergence theorem gives 
rt 



J1 



{\a s {z)Af(z) + b s {z) • V/(z) + H s {z) ■ Vf(z)) Ps (dz) ds. (5.101) 

JC 

Proposition 15 . 25l and the weak convergence results p"(/) — > Pt(f), Po(f) — * Po(f) then complete the 
proof. 1 

Proposition 5.28 Suppose / £ 3 C (or f £ $ if the linear growth condition holds), f = bp-lim/ nfc , 

V/ = bp-limV/nj,, = bp-lim9i9,-/n fc > V/ = bp-lim V/„ fc and R{fn k ) «S fmf/i rifc > lj 

whenever f £ $ c , then for all t 0, 



M t (/„J M t (/) asfc 



(5.102) 



Proof. Set F = sup |||/„J||. Then for every (p,v) £ Cb[0,oo) we have \fn k {z)\ ^ F, pt{F) < oo and 
fc 

hence by the dominated convergence theorem pt{fn k ) — * Pt(/) as /c — > oo. Furthermore I S (F) ^ FA(°\ 
so that I s (f nk ) — > / s (/) an d therefore J Q * J s (.F)ds ^ tFA^°\ A final application of the dominated 

convergence theorem yields J Q I s (f nk )ds — + L I s (f) ds, where we have made use of the continuity 
condition, hypothesis 12.51 to justify the integrability of the LHS with respect to Lebesgue measure. 

Additionally by proposition 15.221 almost surely (p, v) £ Ce[0, oo) n V*, therefore we have for s £ [0, t] 
\\a s {z)Af nk {z) + b s (z) ■ Vf nk (z) + H s (z) ■ Vf nk (z)\ < (io^ + b^ + l)FH% f) ir m (z) (5.103) 



when / £ $ c . We use the bound 



\\a s {z)Af nk {z) + b s (z) ■ Vf nk (z) + H s (z) ■ V/„ fc (z)| < (A<£ + ^ + l)J , fl- <so 7r I11 («) 



(5.104) 



when the linear growth condition holds (recall remark |2.10|I . To deal with both conditions at once set 
Hoo = ^R(f) w ben / £ $ c . Therefore p s ((50^ + &oo + l)FH 00 7r m ) < 00 thus the dominated convergence 
theorem implies that 

p s aa s Af nk + b s ■ Vf nk + H s ■ V/„J 
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p s {ha s Af + b s ■ V/ + H s ■ V/) (5.105) 



i 2 

'^2 



< (kL + boo + IjFH^ sup p s (^ m ) (5.106) 



which is integrable on [0,t] when / G 5 C . and hence 

/ p s (±a s Af nk +b s -Vf nk +H s -Vf nk )ds^ f Ps (}a s Af + b s -Vf + H s -Vf)ds. (5.107) 
Jo Jo 

The integrability of the LHS comes from the continuity of p s (^a s Af rik + b s ■ V/, H . + H s ■ V/„ fc ). 

For the interaction kernel terms with ^= (1, 1) the dominated convergence theorem implies that 

/ K Stijj (Zl, . . .,Zi,dWl, . . .,dWj) < ^2 fn k ( w m) ~ ^ f n k( z n) \ 
J & lm=l n=l J 

-> / K s ,j tj (zi, ...,Zj, dwi, dwj) I f(w m ) - f(z n ) > (5.108) 

lm=l n=l J 

the LHS being bounded by (i + j)K 00 FTr m (zi) ■ ■ ■ 7r m (zi). In consequence a further application of the 
dominated convergence theorem implies that 



/ K s ,i,j(zi, ■ ■ .,Zi,dwi, . . .,dwj) < ^2 fn k (w m ) - ^ fn k (z n ) \ 

^ 3 \7TL — 1 n— 1 ) 



p s (dzi) • • -_p s (dzi) 



J K Sii j(zi, . . .,z i: dwi,. . -,dwj) | ^2 f( w m) - ^2 f P*( dz i) ' ' -Ps(dzi). (5.109) 



We deduce that 
ft 



/ / K Siii j(zi, ■ ■ .,Zi,dwi, . . .,dwj) < fn k {w m ) - ^ fn k (z n ) \ p s (dzi) • ■ ■p s (dz l ) ds 

/ / Ks,i,j(zi, ■ ■ ■ , z l: dw 1: . . . ,dwj) < V" /(w,n) - y^ /(z«) f Ps(dzi) ■ • -^(dz^ds (5.110) 

J J& + > U=l n=l J 

since the integrand on the LHS is continuous by hypothesis 12.81 (and assuming that p s is continuous 
which is almost surely true) and is bounded in modulus by K^F^i + j) sup p s {i^m) 1 which in turn is 

integrable on [0, t]. 

For the self-interaction kernel term i = j = 1, we have for all z £ € by the dominated convergence 
theorem, 

f K sX1 (z,dw){f nk (w) - f nk (z)} ^ f K sX1 (z,dw){f(w)~f(z)} (5.111) 

which is bounded by 2FK^^ if / S $ c and bounded by 2FK oq tt 1u (z) when the linear growth condition 
holds. It follows by the dominated convergence theorem that 



K Sil<1 (z,dw){f nk (w)-f nk (z)}p s (dz)^ / K S:1>1 (z,dw){f(w)-f(z)} Ps (dz). (5.112) 
Under either the linear growth condition or / € 3c the RHS is bounded and integrable in time, so that 



K Sthl {z,dw){f nk {w)-f nk {z)}p s (dz)ds^ / / K StlA {z,dw){f{w)-f{z)}p a {dz)ds. (5.113) 
o Jc 2 Jo J<t 2 

Putting these results together proves the proposition. 1 
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Proposition 5.29 Almost surely, for all f £ §c (or for all f £ $ if the linear growth condition holds) 
and for all t ^ 0, M t (f) = 0. 

Proof. For each t € Q + and neNwe have by proposition 15. 2 71 proposition 15. 121 and Fatou's lemma 

E\M t (f n )\ < liminf E\M^(^ fn )\ = 0, (5.114) 

k 

so that M t (fn) — almost surely. Observe that 

P({M t (/„) ^ for some n and i S Q+}) = (5.115) 

and since every f € $ c (or / G $ if the the linear growth condition holds) is the limit in the sense of 
proposition l5 . 281 of a subsequence f nk (with nt ^ 0) we have that M t (f) — for all t € Q + . Recall from 
corollary 15 . 2 II that almost all the sample paths of (P, u) are continuous in time, and hence M t (f) = 
for all / and for all t ^ except on a set of zero probability. 1 

Definition 5.30 Define the random variables 

v? = sup (u*)- 1 ^,*] and Q? = P$. (5.116) 

Theorem 5.31 The sequence of random variables Q k € D Mb ^ [0, oo) converges in distribution to P 

as k — > oo, where D Mb ^ [0, oo) is equipped with the Skorohod metric d wea k derived from d woa k using 
equation <|2.24[l . 

Proof. Observe that if sup — t| ^ iV - ^ then sup — t| ^ 2iV~3. Now equation Ij5.3|l 

implies that 

/•OO 

^(Q^O < / e" T sup d woak (Qf,P t Ar )dT. (5.117) 

JO O^t^T 

Since d wea k(Q A ', P N ) is a non-negative random variable bounded above by | we have that 



E 



sup d weak (Q?,P t N ) = EV4-"( sup \P$(f n )-P t N (fn)\ A 1) (5.118) 



" n=u 



< 1 A E sup 



M^($ / „)-M f JV (ci> / „)+ / ^(0> / J(Pi v , u f)d s 

(5.119) 

< ^ 4-" ( 2E sup |Mf($/J| + 2A^ig(T + 2)^ (/n) ff- (/n) |||/#|||# ' 



+ |P ( sup |uf - t| > J (5.120) 

^ _^L[(Q(T + 2))*JV-* +N-iQ(T + 2)] + l(T+l)ir-iN-i (5.121) 
v 2 — 1 

on using proposition 15 . 1 2f d) . Thus 

E(rf WC ak(Q A ',P Ar )) ^0 as N -> oo, (5.122) 

it follows that d wcei ] l (Q Nk , P Nk ) — > in probability. Now apply corollary 3.3.3 and remark 3.3.4 of 23 

to deduce that Q Nk P. I 
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6. The Discrete Approximation Scheme 

Theorem 15.311 combined with the construction of sample paths to the martingale problem for jump 
processes (e.g., section 4.2 [23] ) provide us with a technique for generating the sample paths of the 
processes Q N . We start with an initial weak approximation p$ € DJl N to Po at time t — 0. Given the 
sample path generated up to and including time t, the procedure is then to compute r _1 = p N (p^,t) 
from equation (|5.14|) and set pf = pf for s G (f, t + r). With probabilities 

At c 2 N p?{<£) r f K t , i , j (z 1 ,...,z i ,e>)v?( P ?,dz 1 , ...,dz i ) , 



^(pT,*)' P N (P?,ty P N ( P ?,ty nji P N (pf,t) 

for i = 1, . . . ,T and j = 1, . . . ,J, determine p^ +T by rules (1), (2), (3) and respectively where 
rule (1) is to select an element z € £ according to the law It /At and let p^ +T = pf + N~ 1 5 z if this 
is in Wl N , otherwise set p^ +T — pf . Rule (2) is to choose an element of £ according to pf /p^(<£) 
and compute its evolution under the diffusion, external and internal drifts: i.e., choose Z according to 
the N(0, Jdi) distribution and let z' = (n nl (z),j('K^(z),at(z)ZcJ f 1 + b t (z)c]^ 2 ), j^(z)) then let p^ +T = 
pf + N^ 1 (S z > — 5 Z ). Rule (3) is to do nothing: p^ +T = pf . Rule is to select zi,...,Zi according 
to the law 

Kt,i,j (zi,...,z i} £ J )t/ f (pf ,dzi, . . .,dzj) ^ ^ 

^t,i,jy^ii ■ ■ ■ , Hi ) l 
then choose w\ , . . . , w, in accordance with the law 



Kt,ij(z\,...,z i ,&)v i (p t ,dzi,...,dzi) 
e with the law 

K t ,i,j(zi, . . .,Zi,dwi, . . .,dwj) 



K tit j(zi, . . . ,ZiX j ) 



(6.3) 



and set pf +T = pf + N I ^w m — I • In terms of computer implementation, the scheme 

\m— 1 n— 1 / 

may be augmented by the use of majorants and their associated fictitious jumps. The computation of 
efficient majorants for these processes will depend on the precise form of the interactions. 




7. Applications to Gas Phase Dynamics 

The approximation scheme presented in section 6 may be readily applied to the problem of the 
simulation of the formation of particle clusters in a reactor. In this section we point out some of 
the limitations of the method and to what extent they limit practical applications. We also point 
out how the rather general interaction kernels may be constructed to describe a selection of common 
phenomonena. 

One of the major limitations of the current work is that hypothesis 12.81 precludes a purely local 
interaction, i.e., where tVq(zi) = . . . = nji(zi) = ^(wi) = . . . = Tr-^(wj), K t .i.j(zi, . . . ,Zi, dwx, . . . , dwj)- 
a.e. Further restrictions on the diffusion, drift and interaction kernels need to be made to yield the 
corresponding existence results. In the literature, e.g., EH El ED the diffusion coefficients are 
usually assumed to be spatially homogeneous. This rules out some important physical applications, in 
particular we are not able to deal with those situations where the temperature of the gas may vary 
according to its position in the reactor vessel. 

Included within the model defined by equation l|l.l|) is that discussed in |25| . In this paper Turing 
showed by analysing the Fourier series of the concentration that a diffusion-reaction system can exhibit 
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catastrophic instability from a homogeneous equilibrium when non-homogeneous perturbations are in- 
troduced. The catastrophic instability refers to the divergence of the local particle concentration in 
finite time and arises when the magnitude of one of the Fourier coefficients becomes infinite. It follows 
from Parseval's theorem that the particle density is not square integrable for all time and a localized 
version of equation (| 1 - 1 1> with quadratic or higher order terms will not in general have a solution in 
the weak/vague sense. For this reason it seems unlikely that solutions to a localized version of the 
model can be found except in a number of special cases. Let us mention a number of approaches to 
this problem. In |26j . the existence and uniqueness of a suitably regular solution to the Smoluchowski 
coagulation equation was assumed, the stochastic particle approximation could then be compared to 
the solution and convergence established. However Norris has shown that uniqueness need not 
hold for general coagulation kernels even in the homogeneous setting. Deaconu & Fournier |17j use a 
spatially inhomogeneous mass-flow technique to generate solutions to a mollified version of the pure 
coagulation equation. In this case the spatial independence of the diffusion coefficients (and absence 
of drift term) in the domain K rfl allows them to exploit the properties of the heat kernel to make an 
a priori estimate of the local particle density. This estimate is uniform in the quantity e determining 
the mollification, and allows the e — > limit to be taken. The theorem is proved under conditions of 
uniform ellipticity on the diffusion terms. In the discrete case the coagulation kernel obeys a bound 
of the form K(m,m') ^ C(m + to'). Laurengot & Mischler have used weak compactness and a 
diagonal argument together with the results of [2U to deduce an existence theorem under more general 
conditions, in particular fragmentation is taken into account and the uniform ellipticity condition is 
dropped. 

In those cases where no existence result is known for the localized version of equation we can 
make the following heuristic arguments from order of magnitude considerations. Supposing a well- 
behaved limit exists it might be expected that our model approximates the localized version well on 
condition that both the effective range of the particle interaction is large compared to the hydrodynamic 
scale (~ N~ 1/dl ) and tends to zero, i.e., the interaction localizes, as N — > oo. In this way the interaction 
terms correctly detect the limiting particle density as N — ► oo. From the point of view of the simulation 
of molecular dynamics it would be relevant to mollify the interactions on any macroscopic scale over 
which the limiting particle concentration has a suitably small variation, indeed the pointlike particle 
limit assumed in the derivation of equation (| 1 - 1|> is itself unphysical, and would need to be modified in a 
way that prevents infinite particle densities from arising if the molecular dynamics were to be correctly 
described. 

Now let us turn to the possible applications of the model. The individual particles can carry 
additional information that evolve with time and during interactions with other particles; these might 
include discretely changing quantities such as electric charge or the number of active sites, such as in the 
theory of soot formation arising from combustion theory, where free radicals on the surface of the soot 
particles lead to sites where binding to other particles can occur. The number of these sites typically 
decay with time. To incorporate this decay process into the model we include an internal coordinate 
taking values on a closed finite interval which contains all non-negative integers less than the largest 
possible number of active sites for a particle of the given size. The coagulation kernel will depend on 
the number of active sites of each of the potential coagulation partners and the self-interaction kernel 
will describe, possibly amongst other things, the rate at which the active sites decay. 

Often it useful to have additional information about the morphology of the particles under investiga- 
tion. For instance the distribution of surface area is a valuable piece of information for certain chemical 
engineering applications. To account for this we take the particle's surface area as one of the internal 
coordinates and allow this to decay in time according to a size-dependent internal drift term. In this 
way sintering can be quantitatively described in the model. 

The generality of the diffusion and drift terms in equation (jl . 1|> allows us to incorporate the effect 
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of a thcrmophoretic force, this is a force arising from the differential change in pressure acting on the 
particles due to a temperature gradient in the reactor. The thermophoretic force manifests itself as a 
mass-dependent drift term in the modelling equation Hl.lfl . 

Chemical reactions may be modelled by labelling each of a finite number of chemical species by an 
integer, then, provided the chemical reactions proceeds by the law of mass action (if not it may be 
possible to write the reaction as a number of more elementary reactions that do obey the law of mass 
action) a suitable coagulation or fragmentation kernel can be used to implement the reactions. 



8. Conclusions 

In this paper we have constructed a purely stochastic scheme to solve a complicated partial dif- 
ferential integral equation. After imposing a number of hypotheses on the growth of the terms in the 
equation we were able to show the convergence, in a suitable sense, of this scheme to a solution of the 
equation. Our approach differs from many results in the mathematical literature in that we incorporate 
many possible physical processes rather than concentrating on just one or two. In addition we chose 
our hypotheses in such a way so as to maximize the physical usefulness of the resulting scheme; for this 
reason we are satisfied with using mollified interaction kernels if it means, for example, being able to 
model non-trivial diffusion, convection and source terms and realistic reactor geometries. Very few of 
the existing mathematical results allow for temporal inhomegeneity in the parameters describing the 
system. In contrast our result allows for such inhomogeneity, which may typically result from time 
dependent thermal variations in the reactor vessel. 

The aim of this work has been to develop an approximation scheme and associated convergence 
result with wide applicability for the numerical study of molecular dynamics. In the scientific modelling 
literature it is relatively easy to write down a reasonable Monte Carlo scheme depending on many 
variables, though in those cases it is highly unlikely that any existing convergence theorem will apply 
directly. By proving the convergence for a sizable class of such models we have been able to put some 
of these schemes on a firm mathematical foundation. 

One particular objective of this current work was to avoid the necessity of using a hybrid approach of 
Monte Carlo and finite clement methods. A strength of the Monte Carlo approach is its superior scaling 
properties with the dimension of the configuration space compared to finite element methods. It is useful 
therefore to develop a purely stochastic approach to solving these problems. The author believes that 
there is considerable merit in having a number of different approaches to tackling such problems and 
the development of new schemes of all types is welcome, and all the more if the convergence of these 
schemes can be demonstrated. 
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Proposition A.l The space (5c, ||| • |||) defined in definition ^. 1 1\ is separable. 

Proof. Let {p n }^Li be the set of all polynomials on € with rational coefficients and suppose that 
G e C°°(R) is a fixed smooth map satisfying 1(_ <x>,o] ^ ^ ^ oo,i]- ^t 

{9k}f=i = {z~G{\\z\\-r) Vn {z) : n,reN}, (A.l) 

and notice that by using proposition A. 7.1 [23] and the ||| • |||-density of {p n )^Li in the space of all 
polynomials on each compact set of £ that {gfc}^! is a countable set of compactly supported functions 
that is dense in (C£(<£), ||| ■ |||). For each n,m G N, let C„, m = {/ £ i? c : |||/ - .g m || < n^ 1 }, and for 
every pair of integers (n, m) such that C„ !m 7^ choose f n ,m G Cn,m, it follows that 

{/ n ,m : C„, m ^ 0, n, m G N} (A.2) 

is a countable collection of functions that are dense in (5 C , III • III)- 1 



Appendix B. Whitney Extension Result 



Proposition B.l Suppose x G <9f2 and y,z G SI swc/i i/ia£ ||y — ||z — x\\ < HVVwH^ 1 A 5 where 5 
is such that for all w with \\x — w\\ < S we have that \[Vu>(w)\\ > 1. Let 

p = tf-2||*-y||Vw(i/)/||Vw(i/)|| (B.l) 

then the line straight line segments from y to p and from p to z are contained within Q. 

Proof. For 9 G [0, 1] we use Taylor's theorem to compute 

u(y(l-d)+ P 0) < ^(y)-20||z-y||||V^( 2/ )||+20 2 || 2/ -z|| 2 ||VV^|| oc (B.2) 
< -29\\z - y\\(l - \\y-z\\ HVVwIloc) <0. (B.3) 

and 

cj(z{1 - 9) +p9) < uj(z) + 9[{z - y) ■ Vlo( V ) -2\\z- y\\ \\Vu>{y)\\] + 30\\y - z\\ 2 \\ Ww||oo 

+ §0 2 || 2/ -z|| 2 ||VV W || oo (B.4) 

< ^Wy-zUl-fWy-zWWVVuW^^O (B.5) 
since \\y — z\\ ^ \\y — x\\ + \\x — z\\ < ^|| Wt^H^ 1 . This proves the proposition. I 

Proposition B.2 (a) If f G 3c U{/o} then there exists a function f# G C 2 (£#) such that the restric- 
tion satisfies = / and \\\f&\\\& is finite. 

0>) If f G 3c ^ien /# con 6e chosen such that max7r m (supp(/#)) = R(f). 

Proof. If / = /o — 1 then let ft — 1 and the proposition follows. Assume that / G 3c- The result will 
follow from the Whitney extension theorem (|2E], [2H appendix A, [23| proposition A. 6.1). For each 
to G N let (C m = fix T m , €* = co(O) x T m and let F m G C 2 (£ m ) be defined by F m {-) = /(to, •). If 
F m = then let the extension F& G C 2 (£#) be identically zero. In this way part (b) is satisfied. We 
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note that |||/#|||# = sup{||/#(z)||, \\Vf#(z)\\, ||V/#(z)||, \\VVf#(z)\\ : z e (£#, ir m (z) < R(J)} is finite 
by the compact support condition and continuity. Suppose that F m ^ 0, then for £ TO construct 



R°(u,v) = 
R 2 (u,v) = 



F m (u) - F m (v); 

F m {u) - F m (v) - F^(v)(u - v) 







if u 7^ v; 
if u = v, 



F m {u) - F m (v) - F^(v)(u -v)- F£(v)(u -v)(u- v) 



\\u — v\\ 







if u 7^ v; 
if u = v. 



(B.6) 
(B.7) 

(B.8) 



To ascertain the existence of an extension F# € C 2 (£#) it suffices to show that each of these functions 
is (uniformly) continuous on the compact set £ m x £ m . The function R° is clearly continuous, as are 
R 1 and R 2 if u ^ v. Examine continuity at (z,z). Write ir-^ : € m — > Q for the projection onto Q. 
If x = Tifi(z) <Z" dil then there exists a ball, B of radius r > centred at x contained within ft. Let 
C = B x T m then C is a convex open set (in the topology on £ m ) and it follows that if u, v € C with 
m/d then 



\R\u,v)\ 



\f^((1 - t)v + tu) F' m {v)}^-^ At 



<w{F' m ,\\u-v\\)^0 



(B.9) 



as (u, v) —> (z, z) and where the modulus of continuity is defined by equation Ij2.25(l using the operator 
norm. Similarly 



\R 2 (u,v)\ = 



\f^((1 - t)v + tu) F^v)} ^ U) J U V ] (1 - t) dt 



^±w(F£,\\u-v\\)^Q. 



\u — v\\ \\u — v\\ 

(B.10) 

If x £ dfl then there may not exist a convex neighbourhood of z contained in £ m in which case we 
make use of the path given by proposition lB.il Write u = (y,Y) and v = (u>, W), with y,w S f2 and 
Y, W € r m , set p — y — 2\\y — w\\ Vw(y)/|| Vw(y)|| and £ = (p 7 Y) then we have for m/«, 



\R\u,v)\ 



\F^{{l-t)i + tu)-F^{v)]^\ + [F' m {(l-t)v + ti)-F^{v)] {i ] 



\\u — v\ 

^ 5w(F^ n , 3\\u — v\\) — *■ as (it, u) — ► (z, 2) 



||u-«|| 



df 
(B.ll) 



since \\u — £| ^ 2||tt — v\\ and ||£ — u|| ^ 3||u — Similarly 
\R 2 (u,v)\ = 



r (i-t)K((i-i)e +(tt )-f:wi ( "" ( """ ! 



w — u|| \\u — u|| 



+ K((i-^ + ^)-^(«)] 



||u — u|j ||u — u ]| 



di 



(B.12) 
(B.13) 



^(^,311^-^11)^0 as («,«)-» (z,z). 

It follows by the Whitney extension theorem that there exists F# g C 2 (£^) such that -F^le = 
Define f#(m,z) — F#(z) to complete the proof. 
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Appendix C. Particle Selection Measures 



Proposition C.l If P = £^ =1 <5 Zi then f n (P, B\ x- • • x_B„) is f/ie number of injections a : {1, . . . ,n} — > 
{1, . . . , ./V} smc/i i/iai (z Q (i), . . . , z Q (n)) € Bi x ■ ■ ■ x P n . 

Proof. A map [3 : {1, . . . , n+ 1} — > {1, . . . , N} with € B r for r = 1, . . . , n+ 1 and a = n i an 

injection, is either an injection itself, or there exists a unique j € {1, . . . , n} such that a(J) = (3{n + 1) 
with z Q (j) € Bj n B n+ i. Accordingly 

n 

v n+ i(P, B x x • • • x B„ +1 ) = i/ n (P, Si x • • • x B n )P(B n+1 )-J2 M p , Pi x • • • x (S J nB n+1 ) x • • • x B n ) 

3=1 

(C.l) 

with i/i(P, Si) = P(Bi). We show that equation l|5.6[) satisfies this equation. Clearly equality holds 
for v\. We regard 5 rl as the stabilizer of n + 1 in 5„+i. Note that for any a € S n +i either cr € S„ or 
a = (n+ l,j)cr' with a' G 5 n and j = er(n + 1) 6 {1, . . . , n}. Thus 

E C(<0 II p (4n...ns it ) = ^((W+i) n P(B h n...nB ik ) 

o~£S n .+i disjoint cycles w disjoint cycles 

(li,.. M ife) of cr n+l^(zi,...,2 fc ) of cr 

in S n + i in S„ 

n 

+ E E Wn+ 1 nB j )n^(i)n...nB ffl » t)1 ) J] ?(B, 1 n...nB ! J(C2) 

^ — ct'gS^ disjoint cycles 

ff=(n+l,j)o-' jg{i l7 ...,i k ) of o-' 

in S n 

n 

= v n {P, B 1 x ■ ■ • x B n )P(B n+1 ) - E "nOP, Si x ■ ■ • x (Bj n B n+1 ) x • ■ ■ x P n ) (C.3) 

3=1 

= u n+1 (P,B 1 x---xB n+1 ). (C.4) 
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